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Abstract The solar chemical composition is an important ingredient in our under- 
standing of the formation, structure and evolution of both the Sun and our solar system. 
Furthermore, it is an essential reference standard against which the elemental contents 
of other astronomical objects are compared. In this review we evaluate the current 
understanding of the solar photospheric composition. In particular, we present a re- 
determination of the abundances of nearly all available elements, using a realistic new 
3-dimensional (3D), time-dependent hydrodynamical model of the solar atmosphere. We 
have carefully considered the atomic input data and selection of spectral lines, and ac- 
counted for departures from LTE whenever possible. The end result is a comprehensive 
and homogeneous compilation of the solar elemental abundances. Particularly notewor- 
thy findings are significantly lower abundances of carbon, nitrogen, oxygen and neon 
compared with the widely-used values of a decade ago. The new solar chemical compo- 
sition is supported by a high degree of internal consistency between available abundance 
indicators, and by agreement with values obtained in the solar neighborhood and from 
the most pristine meteorites. There is, however, a stark conflict with standard models 
of the solar interior according to hclioseismology, a discrepancy that has yet to find a 
satisfactory resolution. 
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1 INTRODUCTION 



The solar chemical composition is a fundamental yardstick in astronomy, to which 
the elemental abundances of essentially all cosmic objects, be they planets, stars, 
nebulae or galaxies, are anchored. The importance of having accurate solar el- 
er nental a b unda n ces can thus not be overs tated. From the pioneeri i ig effo rts 
of lRusselll (|l929l ). ISuess k UrevI (llQSd'l and IColdberg. Miiller All^ (Il960l^ t o 
the more recen t works oflAnders &: Grevesse ( 19891 ') . Grevesse Sauval ( 1998 1. 
LoddersI (l2003l'l.lAsDlund. Grevesse &: SauvaJ (|2005l ). lGrevesse. Asplund &: SauvaJ 
(|2007l l arid Eodders. Palme k Gaill (|2009l V compilations of the solar system abun- 
dances have found extremely wide-ranging use in astronomy and cosmology. 

There are two independent and complementary ways of determining the solar 
system abundances, each with its pros and cons. Through mass spectroscopy 
of meteorites in terrestrial laboratories, it is possible to directly measure the 
abundance of almost every element and isotope with remarkable precision. The 
most pristine meteorites are the so-called CI chondrites, which have been modified 
least by various physical and chemical processes over the past 4.56 Gyr. Even in 
these meteorites volatile elements have been depleted to various degrees, including 
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the six most abundant elements: hydrogen, hehum, carbon, nitrogen, oxygen and 
neon. As a consequence, it is not possible to rely on meteorites to determine 
the primordial solar system abundances for such elements. This also implies 
that one must measure all meteoritic abundances relative to an element other 
than hydrogen, traditionally chosen to be silicon. Prior knowledge of the solar 
photospheric Si abundance is therefore required in order to place the m eteoritic 
abundances on the same absolute scale as the Sun (jSuess fc Urevlll95fil l. 

With the exception of a general ~ 10% modification due to diffusion and grav- 
itational settling (Sect. 13. lip and depletion of lithium and possibly beryllium, 
today's photospheric abundances are believed to reflect those at the birth of the 
solar system. They are of profound importance for the understanding of our own 
solar system and the Sun's interior structure and evolution. Unfortunately, pho- 
tospheric abundances cannot be determined with the same accuracy as those from 
meteorites. In particular, very little isotopic abundance information is available. 
Furthermore, photospheric abundances are not observed but inferred. Thus the 
solar spectrum must be interpreted using realistic models of the solar atmosphere 
and spectrum formation process. In addition, the required atomic and molecular 
data are rarely precise enough to rival the accuracy of meteoritic measurements. 
In spite of the remaining uncertainties in solar photospheric determinations, hav- 
ing access to the most accurate solar abundances possible is crucial, not only for 
understanding our own Sun but also in order to make meaningful comparisons 
with other stars (for which meteoritic information is obviously unavailable). The 
Sun is also the only possible reference for the interstellar medium, Hii regions 
and other galaxies when studying the elements depleted in meteorites. 

In the decade since the a ppearance of the w i dely-u sed compilation of the so- 
lar chemical composition by lGrevesse &: Sauva l (ll998l ). the abun dances of a large 

numb er of elements have been revised, several rather severely (e.g. Asplund. Grevesse &: Sauval 
iooi). This is true in particular for C, N, O and Ne. There are three main reasons 



for the changes: improvements in atomic and molecular transition probabilities, 
the advent of 3D hydrodynamical solar model atmospheres and relaxation of 
the assumption of local thermodynamic equilibrium in the spectral line forma- 
tion. Several of the proposed abundance alterations have not gone unchallenged, 
however. It is therefore timely to now critically assess the reference solar chemi- 
cal composition and make recommendations regarding the most accurate values 
available to the astronomical community. After recapitulating the necessary in- 
gredients in any solar abundance analysis (Sect. [2]), we discuss in detail the solar 
chemical composition on an element-by-element basis (Sect. [3]). We also present 
a comparison between photospheric solar abundances and those inferred by alter- 
native means, such as via meteorites and stars or gas in the solar neighborhood 
(Sect. H]). We conclude by summarizing the main points and outlining a few open 
issues for the future (Sect. [5]). 



2 INGREDIENTS FOR SOLAR ABUNDANCE ANALYSIS 
2.1 Observations 

Analyses of the solar photospheric chemical composition make use of the ele- 
mental fingerprints present in the solar spectrum in the form of spectral ab- 
sorption lines. There are several solar atlases of both intensity and flux spectra 
in the ultraviolet, optical and infrared regions. The most commonly used opti- 
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cal disk-center intensity spectra of the quiet Sun are the so-called Jungfraujoch 
or Liege atlas (after its locations of observation and productio n, respectively; 
Delbouille. Roland &: Neven 19731 ) and the Kitt Peak solar atlas ( Neckel &: Labd 
1984l ;i. The spectral resolving power of the latter is slightly higher while the for- 
mer is less affected by telluric absorption due to the higher observing altitude. 
There are a,lso correspondiiig infr ared disk-center intensity atlases observed from 
Kitt Peak ( Delbouille et al. 198ll) and from the ATMOS ex periment flown on the 
space shuttle ( Abrams et al. 19961 . Farmer &: Norton 19891 ). Finally, the solar flux 
spectr um has recently been re-reduced by Kurucz ( 20061 ): see also Kurucz et al 



(|l984l ) "or an earlier rendition of the same Kitt Peak solar data set. All solar 
atlases agree very well with each other except for spectral regions afflicted by 
telluric features. Thus the quality of observations is in general not a source of 
significant error in solar abundance analyses. 



2.2 Atomic and molecular data 

A critical ingredient in any solar abundance analysis is of course the atomic 
and molecular line input data. Of these the transition probability is the most 
obvious, but there are a raft of other required components that are not al- 
ways as well-determined as one would hope. These include line broadening, 
hyperfine and isotopic splitting, dissociation energies and partition functions. 
Much laudable work h as been done in improv i ng the q /-val u es for a large num- 
ber of transitions (e.g. Johansson et al.l l2003l . iKuruca Il992l . iLawler et al.l I2OO6I . 
Sobeck. Lawler SnedeiJ 200?! ) but still several elements and lines have surpris- 
ingly uncertain transition probabilities. Another important advancement in re- 
cent years is that it is now possible to accurately com pute the self-broadening 



(also known as van der Waals broaderiing) of metal lines (lAnstee fc 0'Maralll995l . 
Barklem &: Aspelund-Johanssonll2005l . Barklem. Piskunov &: O'Marall2000al ) , which 
has made the classical Unsold broadening and its ad-hoc enhancement factors 
largely obsolete, at least for neutral species. 

The need for atomic data increases dramatically when relaxing the LTE as- 
sumption. In non-LTE calculations one needs transition data not only for the 
line in question but in principle for all other lines of the element. In addition, 
photo-ionization cross-sections play a crucial role. There have been notewor- 
thy impr ovements in this res pect in recent years through the Opacity and Iron 
Projects (|Badnell et al.l 120051 ) for example, but unfortunately the lack of quan- 
tum mechanical calculations for most elements preclude detailed non-LTE in- 
vestigations. The most unsatisfactory situation, however, relates to collisional 
cross-sections for excitation and ionization, for which laboratory or theoretical 
calculations are completely missing bar for a few elements and transitions rele- 
vant for spectroscopy of late-type stars (e.g. Barklenj 2007al . Belyaev fc Barklem 
2003l'l. In the absence of re alistic ev a luatio ns, classical recipes such as those of 
van Regemorter (|l962l ) and IPrawinl (|l968l ) are frequently used but cannot be 
expected to be more than order-of-magnitude estimates at best. Often scaling 
factors to these formulae are used, sometimes calibrated by achieving an opti- 
mal overall fit to the various spectral features. It should remembered that such 
a procedure may reveal less about the true collisional cross-sections than about 
other shortconi ings in the modelling of the stellar atmospheres and line formation 
(|Asplundll2005l l. 
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2.3 Solar atmospheric and line formation modelling 



Extracting information about the Sun's chemical composition from the solar spec- 
trum requires realistic models of the solar atmosphere and the line formation pro- 
cess. Traditional solar and stellar abundance analyses employ one-dimensional 
(ID) model atmospheres assuming time-independence and hydrostatic equilib- 
rium. Such models come in two flavours: theoretical and semi-empirical. The 
former are constructed assuming a constant total (radiative plus convective) en- 
ergy flux through the atmosphere, which dictates the temperature stratification 
given the adopted equation-of-state and opacities. Convective energy tra nsport 
is normally described by the mixing length theory ( Bohm-Vitense 19581 ). The 
radiative flux is determined by solving the radiative transfer equation, typically 
under the simplification of local thermodynamic equilibrium (LTE). To achieve 
a realistic temperature structure, treatment of the combined effects of spectral 
lines ('line blanketing') must be as complete as possible. Th e most widely- used 
theoreti cal ID model atinosphe res for solar-type stars are the lKurucj(jl99.ll ) and 
MARCS (jCustafsson et ahllioosl ) grids of models. 

Semi-empirical models are often the preferred choice for solar spectroscopists. 
In these models, temperatures are inferred from depth-dependent observations 
such as the center-to- limb variation of the continuum and/or lines with differ- 
ing formation heights. The resulting atmospheric structure is thus sensitive to 
the ingredients involved in the inversion process, such as LTE for the contin- 
uum and line formation, and the adopted line data. Hydrostatic equilibrium 
is still assumed in semi-empirical models, but the flux constancy criterion does 
not need to be obeyed nor is it necessary to estimate convective energy trans- 
port. Several sein i -empi rical solar models are on the market, inc luding the 
Holweger fc Miilled (Il974l') val3c (jVernazza. Avrett k Loeserl [l976l ) and MISS 
( Allende Prieto e t al.^ ^2001 ) models. The favoured model of solar abundance 
aficionados is the Holw eger &: Miilleii model, partially because its temperature 
stratification is expected to be an accurate representation of the real photosphere 
- but also partially just by h abit. Indeed it is soin ewhat surprising that the val3c 
model or its successors (e.g. iFontenla et "al]l2006l ) have not caught on more within 
the solar abundance community, since th ey are built on simila r premises but with- 
out the LTE restriction. The venerable Holweger &: Miilled model in fact dates 
back more than forty years, since t he 19 74 version was simply a new pressure- 
integration of the original lIolwege3 ( 19671 ) model using updated equation-of-state 
and continuous opacities. It was constructed to match the observed continuum 
center-to-limb variation as well as the line depths of some 900 lines of varying 
strengths, under the assumpti on of LTE . Due to the limited spectral resolution 
of the solar spectrum used by iHolweged . one would expect the temperatures to 
be overestimated by about 50 K in the line-forming region. Interestingly, the cor- 
rected version of such a temperature structure wou ld be more similar to the miss 
semi-empirical model (j Allende Prieto et al.l [20011 ): both semi-empirical models 
have a shallower temperature gradient than ID theoretical model atmospheres, 
as seen in Fig. [H It should be borne in mind that there exists an ambiguity when 
employing any ID model: while many simply adopt temperatures, pressures and 
densities directly from the published model, a more correct procedure is to per- 
form a pressure-integration using only the temperature structure as a function of 
continuum optical depth, in order to be consistent with the opacities computed 
with the spectrum synthesis code one employs. 
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More recently it has become tractable to perform 3D, tin ie-d epende nt, hydro- 
dynamical simulations of the s t ellar surface convection fe.g. lAsp lund et al. 199S , 
Frevtag. Steffen k Dorchll200l IStein Nordlundll 19981 . IVogler. Bruls k Schiissler 



Here the standard conservation equations of mass, momentum and energy 



are solved together with the 3D radiative transfer equation in a small but rep- 
resentative volume of the stellar atmosphere. Radiative energy exchange with 
the gas in the surface layers is crucial, because this is what drives the convec- 
tive motions. Because of the much more computationally demanding nature of 
such hydrodynamical modelling, approximations must be made in the radiative 
transfer treatment. Rather than solving for many thousands of frequency points 
as regularly done in theoretical ID hydrostatic models, the wavelength varia- 
tion of opacity is represented by a small number of typical opacity bins (4-20) 
in order t o compute the radiative heating/cooling rate that enters the energy 
equation ( Nordlund 19821 ). Detailed comparisons with the monochromatic so- 
lu tion reveal that this scheme is sur prisingly accurate. The reader is referred 
to iNordlund. Stein fc Asplundl (|2009l 'l for a recent review of the physics of solar 
surface convection and the numerical details of its simulation. 

A few different 3D hydrodynamical solar model atmospheres have been con- 
structed and employed for abundance purposes (e.g. I Asplund et al.ll2000bl .l ICaffau et al 



2008al . iTrampedach et al.l [2009'). While these are computed with different codes, 



they are based on many of the same physical assumptions and approximations. 
It is therefore not too surprising that the resulting atmospheric structures are 
similar. The 3D models have been computed using realistic and comprehensive 
equation-of-state, continuous and line opacities. The numerical resolution is suf- 
ficiently high to adequately resolve the granulation, albeit obviously not high 
enough to reach down to the microscopic viscous dissipation length-scales. The 
reason that these so-called 'large-eddy' simulations are nevertheless quite success- 
ful in reproducing the real photospher e is that the mean stratification is largely 



indep endent of numerical resolution ( Asplund et al. 2000al . Stein &: Nordlund 



19981 ). Spectral line formation is heavily biased to the upflows, which have a 



low degree of turbulence in comparison to the downflows, due to the divergent 
flow. Obviously none of the mixing length parameters required in ID models 
enter into 3D simulations, as the convective energy flux is instead a natural con- 
sequence of the hydrodynamics. 

Having a suitably realistic solar model atmosphere, the next step when trying 
to infer the solar chemical composition is to model how the spectrum is formed. In 
order to derive an abundance of an element, solar spectroscopists either measure 
the total observed line strength as defined by the equivalent width, or directly 
compare the observed and theoretical spectra through the fitting of line profiles. 
Finally, to compute the line strength requires knowledge of the level populations. 
In LTE, these follow straightforwardly from the local temperature and electron 
pressure using the Boltzmann and Saha distributions. In general this rather 
severe simplification cannot be expected to be valid, and the non-LTE case must 
instead be considered. While non-LTE is a general term covering all situations 
when LTE is not valid, in practice all calculations referred to herein assume 
statistical equilibrium, i.e. the level populations do not change with time. Non- 
LTE thus requires a simultaneous solution of the rate equations for all relevant 
levels and species togeth er with the rad iative transfer equation for all relevant 
wavelengths. We refer to Asplund ( 20051 ) for a detailed account of the principles 
of non-LTE and the most important non-LTE effects encountered for late-type 
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stars like the Sun. 



2.4 Observational constraints on solar modelling 



Since more sophisticated modelling does not automatically imply a more realistic 
outcome, it is of paramount importance to carefully test the predictions of the 
3D model atmospheres against an arsenal of observational diagnostics. This is 
especially true given that ID model atmospheres more often than not fail the very 
same tests. Current generations of 3D models are very successful in reproducing 
the observed solar granulation topology, typical len gth- and time-scales, convec- 
tive velocities and intensity b rightness contrast (e.g. iNordlund. Stein &: Asplund 
20091 . IStein fc NordhmdHToOsi ^. clearly none of which ID models are able to pre- 
dict. 

Another important constraint is the center-to-limb variatioii of the contin- 
uum as a function of wavelength. Avres. Plvmate &: Kelleil ( 20061 ) have criticized 
the 3D hydrodynamical model of lAsplund et al.l (|2000bi r oft used in recent so- 
lar abundance analyses, for having too st eep a temperature gradient as judged 
by th e observed center-to-limb variation. iKoesterke. Allende Prieto &: Lambert 
( 20081 ) reach a somewhat milder conclusion, presumably because they performed 
spectral synthesis with the full 3D model rather than with only the spatially- 
av eraged version used b v lAvres. Plvmate Kelleii . The 3D coSbold solar model 
of ICaffau et al.l (j2008al ) appears to perform better in this respect (H. Ludwig, 
private communica tion). In the present w' ork, we make extensive use of the 
3D solar model by Trampedach et aP ( 20091 ). which has been constructed with 
an improved treatment of the radiative transfer and updated opacities. Be- 
ca use of its some what shallower temperature stratification than the 3D model 
of Asplund et al.l (Fig. [1]), the new 3D model satisfies the center-to- limb vari- 
ation constraint very well, as illust rated in Fig. [2l Indeed it outperforms all 
ID models, even the semi-empi rical (Holweeer fc Miiller 19741 ) model, which was 
designed to fit this diagnostic ( Pereira et al. 20091 ) . 

The wings of the hydrogen Balmer lines are also sensitive tracers of the temper- 
ature stratification in the deeper layers of the photosphere. In recent years, signif- 
icant improvements have b een made in the treatment of the broadening of H lines, 
especially self-broadening (jAllard et al ] l2008l : lBarklem. Piskunov 'Marall2000bl ). 
It is still unclear whether noticeable non-LTE effects can be expected in the 
solar atmosphere, since availa ble calculations for the cross-sections of inelastic 
H+H collisions are uncertain ( Barklem 2007bl ) . With ID theoretical model at- 
mospheres, the H lines imply too low Tcfj for the Sun by 50 — 100 K, and have 
the wrong line shap es when the most accurate line broadening data are used 



(|Barklem et al.ll2002l V Such problems can be partly cured by playing with the 



mixing length parameter s, however. In contrast, the H lines suggest that the tem- 
perature gradient_ofthejHolweger_&^ ( 1974 ) semi-em pirical model is too 



shallo w. iPereira et al.l (|200S ) found that the 3D solar model of lTrampedach et al.l 
( 20091 ) reproduces the Ha and H/? lines very well in LTE, doing so significantly 
better than the ID models they investigated, but without requiring the tweaking 
of any free parameters. The coSbold model has to our knowledge not yet been 
tested against H lines. 

As seen in Fig. [3l with 3D hydrodynamical models for the solar atmosphere it 
is now possible to achieve highly satisfactory agreement with observed profiles of 
typical weak and intermediate strong lines (the wings of strong lines reveal more 
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about the pressure broadening data than the employed model atmosphere). It is 
important to remember that this is achieved without invoking the free parameters 
necessary in any ID analysis (micro- and macroturbulence) in the line formation 
calculations. Such line broadening can thus be explained as the result of the 
Doppler shift arising from th e convective motions with a. smal ler contribution 
from the solar oscillations (e.g Nordlund. Stein &: Asplundl 2009 . and references 
therein). Indeed, even the observ ed line shifts and asy mmetries are very well 



reproduced with 3D models (e.g. lAsplund et al. 2000^^1. The c oSbold solar 



model produces good overall line profiles (e.g. Caffau et al. 2008al ) but whether 
this extends also to line asymmetries has not yet been investi gated. 



In summary, the 3D hydrodynamical model atmosphere (jTrampedach et al 



20091 1 that we draw results from in this review outperforms ID model atmospheres 



in the available observational tests, from granulation topology and center-to-limb 
variation to H lines and the detailed profiles of metallic lines. For the purposes 
of photospheric abundance determinations, the 3D solar model employed here 
appears to be a very realistic representation of the solar photosphere. No large 
differen ces are expected fro m using our 3D model or the corresponding co5bold 
model (jCaffau et al.ll2008al l given their overall similarities. 



3 PHOTOSPHERIC ABUNDANCES 



Our recommended solar photospheric elemental abundances are listed in Table 
[H which also provides the corresponding meteoritic values for CI chondrites to 
be discussed in Sect. 14.11 [Throughout this review, we adopt the customary 
astronomical scale for logarithmic abundances where hydrogen is defined to be 
log Eh = 12.00, i.e. log ex = log{Nx/Ni{) + 12, where A'x and A^h are the number 
densities of element X and hydrogen, respectively.] Fig. [4] shows how the solar 
abundances vary with atomic number, illustrating several key features of nuclear 
and stellar physics: the high primordial H and He abundances; the fragile nature 
of Li, Be and B; the relatively high abundance of elements involved in stellar H-, 
He- and C-burning, modulated by the odd-even effect and the a-capture effect; 
the high nuclear binding energy and near nuclear statistical equilibrium for the 
Fe-peak; and the production of the heavy elements through successive neutron 
capture wit h peaks around the magic nuclei possessing closed neutron shells (e.g. 
Pagellll997l ). 



In this review, we have attempted to reanalyse the solar abundances of (nearly) 
all elements in a homogeneous manner with the best possible atomic data and 
state-of-the-art solar modelling. As a guiding principle, we have been very dis- 
cerning when selecting the lines for each element, since inclusion of dubious lines 
only increases the abundance scatter and tends to skew the results towards higher 
abundances due to blends. The analysis has been carried out using several ID 
and 3D model atmospheres and with non-LTE effects accounted for whenever 
possible, all done with the same well-tested computer codes. Unless specified 
otherwise, the results prese nted here have been base d on the 3D hydrodynamical 
solar model atmosphere of iTrampedach et al. I (l2009l l. Below we present in some 
detail how the abundances were derived but the full description of the analysis 
(including line lists with all the relevant data for the transitio ns) will appear 



in a fo rth coming series of articles in Astronomy & Astrophysics (| Asplund et al. 



2009al lblR iGrevesse et"^l2009l . ISauval et al]l2009l . IScott et al.ll2009bl ). 
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There are a wide range of potential sources of error in solar and stellar abun- 
dance analyses, from inaccurate input data for the transition (^/-values, line 
broadening etc) to difficulties in estimating line strengths (finite S/N, continuum 
placement, blends etc), to inadequacies in the atmospheric and line- formation 
modelling (ID, LTE, mixing length and microturbulence parameters, continuous 
opacities, etc). No consensus exists in the solar abundance literature on how the 
uncertainties should be quantified. Some authors have used the standard devi- 
ation arising from the chosen sample of lines to represent the total error for a 
given species whilst others have employed the standard error of the mean or some 
different measure altogether. Previous compilations of the solar chemical com- 
position contained a mixture of different error estimates, since the values were 
taken from sources with differing modus operandi. 

We have attempted to quantify three possible systematic errors introduced 
by the modelling: mean atmospheric stratification, atmospheric inhomogeneities 
and departures from LTE. The first uncertainty has been estimated by taking half 
the difference b etween the results w i th th e spatially-averaged 3D model (denoted 



<3D>) and the lHolweger &: Miillerl (jl974l l model. The second has been evaluated 



as half the difference between th e full 3D results and th ose from <3D>. The in- 
clusion of the theoretical ma rcs ( Gustafsson et al. 20081 ) and semi-empirical MISS 



(|Allende Prieto et al.l [2OO1I 1 models serves as a further check on these two sys- 



tematic errors. Due to the general lack of non-LTE calculations for the majority 
of elements, the non-LTE uncertainty is more difficult to quantify. Somewhat 
arbitrarily, in most cases we have chosen half of the predicted non-LTE abun- 
dance correction as an estimate of the error. Occasionally we have employed 
the differen c e betw een the cases with and without inelastic H collisions through 
the Drawin (Il968l ) formula, whilst taking into account the recommendations by 



Asphmdl ^200,^ ). We have used an arbitrary, minimal estimate of the possible 



errors in the non-LTE corrections of 0.03 dex. In some cases this might over- 
estimate the non-LTE systematic error, whilst in others (e.g. Til) it probably 
underestimates the error; there is unfortunately no other consistent way to in- 
clude such unknown errors. Estimating possible systematic errors arising from 
inappropriate line input data has not been attempted here; relative errors in 
^'/-values will at least be accounted for in the dispersion of the lines. When 
the elemental abundance is determined by only one or very few lines, we have 
attempted to estimate an additional error arising from uncertainties due to con- 
tinuum placement and possible blends. The systematic errors have been added 
in quadrature with the statistical error to estimate the total uncertainty, which 
is used throughout the article unless specified otherwise. The statistical error 
has been computed from the weighted standard error of the mean, with weights 
assigned to each line according to uncertainties in the continuum placement and 
known or suspected minor blends. 

3.1 Lithium, Beryllium and Boron 

Lithium: In the Sun, Li is depleted by about a factor of 150 compared with 
the already small meteoritic value. Even the Lii resonance line at 670.8 nm is 
therefore very weak. The line is blended by an assortment of CN and Fei lines, 
which makes the Li abundance determination very challenging. Because the 
transition arises from the ground state of a minority species, the line is sensitive 
to both the temperature and photo-ionizing UV radiation field (| Asolundl [2005I . 
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Kiselmanlll997ll. The recommended valu e Kiven in Table [T] stems from the analy- 
sis of Miiller. Pevtremann &: de La Reza ( 19751 ). which has been c orrected for the 



effects of atmospheric inhomogeneities a nd 3D non-LTE effects (jAsplund et al 



19991 . iBarklem. Belvaev fc Asplundl 120031 ). The solar Li depletion cannot be ex- 
plained by standard models based on the mixing len gth theory of convection but 
requires additional mixing below the convection zone (jBrun. Turck-Chieze &: Zahn 
19991 : ICharbonnel fc Talonlliooil ). 

Beryllium: Early work implied that Be w as depleted in the sola r convection 
zone by about a factor of tw o (ff hmielewski, Brault &: Miiller 19751 ). This inter- 
pretation was challenged by iBalachandran fc Bell ( 19981 ). who a rgued that one 
must take into account the long-suspected missing UV opacity (e.g. lShort &: Hauschildt 
2OO5I . and references therein) in synthesizing the Be 11 resonance lines at 313 nm. 
They estimated the amount of missing opacity by requiring that the nearby OH 
electronic lines would return the same O abundance as the OH vibrational lines 
in the infrared, which increased the derived photospheric Be abundance to the 
meteoritic value. The required extra opacity caii be la rgely attributed to photo- 
ionization of Fe I (jBell. Balachandran &: Bautistal[200ll ). The lBalachandran &: Belli 
analysis was corroborated by Asplund ( 20041 ). who found that the basic conclu- 
sions are independent of the adopted model atmosphere, w hether ID or 3 D. The 
Be II lines are not sensitive to non-LTE effects in the Sun (|AsDlundll2005l ). 

Boron: The problem of the missing UV opacity may also impact the solar 
B abundance determinat ion, which can r ealisti cally only be based on the Bi 
249.7 nm resonance line. Cunha &i SmithI (jl999l ) have re-analysed the line with 
careful consideration of the impo rtant Mgl photo-ionization cross-sections. The 
effects of departures from L TE (IKiselman &: CarlssonI Il996l) an d 3D hydrody- 
namical model atmospheres ( Asplund. Grevesse Sauval 20051 ) are small and 
opposite in sign. There is still a large uncertainty attached to the resulting abun- 
dance, log ee = 2.70±0.20, due to the difficulty in analysing this crowded spectral 
region. Not surprisingly, B does not appear depleted in the solar photosphere. 



3.2 Carbon 

The solar carbon abundance can be inferred from a wealth of indicators, from low- 
excitation forbidden and high-excit ation permitted atomic lines to various molec- 
ular transitions of CH, C2 and CO (Lambert 19781 ). This fortunate circumstance 
has not translated into a well-determined abundance over the years as the difi^er- 
ent diagnostics have often yielded discr epant results. The prelimi nary abundance 
logec = 8.56 ± 0.04 recommended bv lAnders Sz Grevessd (19891) eventually be- 



came 8.60 ±0.05 in the final analysis of I Grevesse et al. ( 1 99 ll )7 Grevesse Sauval 
(119981) based their v alue o f 8.52 it 0.06 on an unpublished analysis employing a 
Holweger fc Miiller (Il974l ) model with a modified temperature structure, in an 
attempt to remove existing trends with excitation potential and equivalent widths 
for C, N, O and Fe. Over the past decade t he solar C abundance has undergone 



an ev en more drastic downward revision. lAllende Prieto. Lambert Sz Asplund 
(|2002l ) reanalysed the [C i] 872. 7 nm line and found t hat the combination of a 
3D hydrodynamical solar model ( Asplund et al. 2000bl ) and an updated 5/- value 
brought the abundan ce down to 8.39 it 0.04. The same value was obtained by 
Asplund et al.l ( 20051 ). who also employed Ci, C2 and CH electronic and vibra- 
tion lines using the same 3D model. Convincingly, all indicators gave consistent 
results, further strengthening the case for the reality of the low C abundance. 
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In sharp contrast, the abundances estimated with the Holweger &: Miiher ( 19741 ) 
model atmosphere disagreed by « 0.2 dex. The presence of temperature inhomo- 
geneities and a cooler me an temperature struc t ure in the 3D model compared with 
the previously employed iHolweger fc Miillei] (|l974l ) model reduced the deduced 
abundance from the molecular lines in particular, while significant departures 
from LTE had a similar effect on the C i lines. Using weak lines of the extremely 
temperatur e -sensi tive CO molecule and a self-consistently derived O abundance, 
Scott et al. corroborated this finding. They also determined an isotopic 

ratio of i^C/i^C = 86.8 ± 3.8. 

The value we recommend h ere, l ogec = 8.43 it 0.05, is slightly higher than in 
Asplund. Grevesse &: Sauval ( 20051 ) an d stems from the mea n of the 3D-based 
results for [Cl], Cl, CH and C2 lines (jAsplund et al.1 l2009al ') . As illustrated in 
Table O the various C indicators imply highly uniform abundances. With the 
exception of a correlation between derived abundance and line strength for Ci, 
there are no trends with transition properties such as excitation potential or 
equivalent wid t h for any of the species. A similar trend was also present in 



Asplund et aP (|2005l l. which was attributed to underestimated non-LTE effects 



for the stronger lines ( Fabbian et al. 20061 ). We note that only ID non-LTE 



abundance corrections have been employed for the Ci results while the depar- 
tures from LTE are expected to be exacerbated in the presence of atmospheric 
inhomogeneities ( Asplund 20051 ): a full 3D non-LTE study for C would clearly 
be worthwhile. Given the importance o f this abundant element, it is essential 
that the results of Asplund et al.l ( 2009al ) are confirmed by an independent group 
using an independen t 3D model, akin to the analysis of the solar O abundance 



(jCaffau et ahlboOSal ). 



3.3 Nitrogen 

The nitrogen abundance can be determined from both atomic and molecular 
lines, although the forbidden [Nl] li nes are too weak in the solar spectrum to 



be measurable (Grevesse et al 



m 



199d) . A preliminary N abundance was given 



Asplund. Grevesse Sz Sauval ( 20051 ) using high excitation Ni l ines and NH 



vibra tion-rotation lines in the infrared while in the final analysis (| Sauval et al 



20091 ) also NH pure rotation lines and a multitude of transitions from various 
bands of CN were considered in the framework of an improved 3D model. The 
non-LTE eff ects for N i ar e relatively small: ~ —0.05 dex when ignoring H col- 
lisions fCaff au et ai]l2009l ). In their own analysis of the se very weak Ni lines 
using their coSbold 3D solar model, ICaffau et al. I (I2OO9I ) obtained a Nl-based 



abundance which is 0.05 dex higher than the corresponding value in lSauval et al 



(20091 largely due to the selection of lines. They included many lines that we 
consider to be of dubious value due to known or suspected blends. This is also 
reflected in the dramatically different standard d eviations for t he N i lines in the 
two studies: a = 0.04 and 0.12 dex, respectively. Sauval et al. ( 20091 ) found very 
good agreement between abundances inferred from Ni and NH lines, with very 
small dispersion in the NH results. The different CN bands also gave very simi- 
lar results when empl oying a self-consiste ntly derived C abundance with the 3D 
model. Unfortunately, Caffau et al. ( 20091 ) did not consider molecular lines and it 
is therefore not known how well their 3D model performs in achieving consistent 
results for atoms and molecules. It is sometimes stated in the literature that the 
molecular lines are more temperature sensitive than atomic lines, but in fact the 
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NH lines are less so than these highly excited Ni lines. 

The N abundance from Sauval et al. ( 20091 ) th at we recommend he r e, log 



7.83 ± 0.05, is 0.09 de x lower than sugges t ed by iGrevesse &: Sauvai ( 19981 ) and 



1 .66 ± U.Ub, IS U.uy de x lower tnan sugges t ed by 
0.22 dex lower than in lAnders fc Grevessd ^Wm . 



3.4 Oxygen 

Oxygen is the most abundant element in the cosmos not produced in the Big 
Bang, and the third most common overall. Oxygen and its isotopes are key 
tracers of the formation and evolution of planets, stars and galaxies and as such 
it is one of the most important elements in all of astronomy. Nevertheless, or 
perhaps consequentially, its abundance often seems to be in dispute with the 
solar O abundance being no exception. The favoured solar content has dropped 
almost precipitousl y over the past t wo decades from logep = 8 . 93 ± 0.04 in 
Anders fc Grevessd (|l989l ) to 8.66 in lAsplund. Grevesse k Sauvall (l2005ll. The 



value we recommend here, logeo = 8.69 ± 0.05, is taken from lAsplund et al. 



(j2009d ). who presented a 3D-based study of the available diagnostics: forbidden 
[O i] and permitted O i transitions as well as vibration-rotation and pure rotation 

li nes of OH. 

lAllende Prieto. Lambert Asplund ( 200ll ) reanalysed the [Oi] 630 nm line 
with a 3D model atmosphere. Based on the line profile they demonstrated that 
the li ne is blended by a Nil line, which was subsequently confirmed experimen- 
tally ( Johansson et al. 20031 ). In combination with the 3D model, recognition of 
the blend led to a substantial low ering of the s olar O abundance compared with 
the canonical value at the time. lAvred (|2008l l revisited the 630 nm line with a 
single snapshot of another 3D model but found a 0.12 dex higher O abundance, 
partly because he allo wed the Ni I q /-v a lue to vary freely in spite of the new accu- 
rate measurement by [Johansson et al. (|2003l ). iGaffau et al.l (|2008al ) used several 
snapshots from the same C05B0LD 3 model but instea d ^ fixed the Ni contribu- 
tion assuming th e Ni abundance from Grevesse Sz Sauyal (Il998l ) t hereby confirm- 
i ng th e results of Allende Prieto. Lambert &: Asplund ( 200ll ) and Asplund et al. 
A more self-consistent approach is t o first determine th e Ni abundance 



before considering the [Oi l line, as done bv IScott et al.l (|2009al ) wi th the same 
3D model as employed by Asplund et al.l ( 20041 ) . As explained by Scott et al. 



it is then possible to very accurately predict the Nil contributio n to the overall 
line strength, independent of the reference atmospheric model. lAsplund et al. 



(l2009c ) updated this analysis and found that the effects of the increased Ni 
proportion and their shallower model tern perature gradient in comparison with 
lAllende Prieto. Lambert &: Asplundl ( 200ll ) compensate for one another, leaving 
the abundance at logeo = 8.66. This ab undance is also consistent with the ce nter- 
to-limb variation of the 630 nm feature ( Pereira. Asplund Sz Kiselman 20091 ). 

In addition to the 630 nm line, there are two other [Ol] transitions that can 
be utilized. The 636 nm line is located in the midst of a very broad Cal auto- 
ionization line as well as blended by CN lines; the latter's contribut ion can how- 



ever be accu rately predicted from o ther neighboring CN lines. Both lCaffau et al 



(|2008al ) and lAsplund et al.l (|2009d l :bund that this line returns a slightly higher 
abundance than the 630 nm line, but still consistent within the uncertainties. 
Melendez Asplundl (|2008l l added the [O l] 557 nm line, whose slightly excited 



lower level makes it largely independent of the model atmosp here. Unfortunately 
this transition is badly blended by C2 lines, whose contribution lMelendez &: Asplund 
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attempted to estimat e using nearby C2 line s with the same Hne strengths and 



excitation potential. lAsplund et alj ( 2009d ) obtained an abundance from the 
557 nm line exactly in between those from the other two forbidden transitions. 
The three [O i] lines imply a solar O abundance of log eo = 8.70 it 0.05 when due 
consideratio n is given to the blends. 

Recently, Centeno &: Socas-Navarrj ( 20081 ) introduced a novel method to es- 
timate the O abundance, namely spectropolarimetry of the [O i] 630 nm line in 
a sunspot. Because of the different polarization of the [Oi] and Nil lines it is 
possible to disentangle the Nil contribution. After consideration of the fraction 
of O tied up in CO in the co ol environme r its of s unspots, they obtained a high O 
abundance of logep = 8.8 6. Scott et al. ( 2009al ) have however pointed out that 
Centeno &: Socas-Navarrj applied an outdated gf -value for the [O i] line. Fur- 



thermore, they showed that when using an improved Ni abundance and an alter- 
native CO treatment one finds a significantly lower O abundance: logeo = 8.71; 
a further reduction is obtained when employing an alternative sunspot model. 

There are a half-dozen or so permitted, high-excitation O i lines useful for so- 
lar abundance purposes. Of these the relatively strong 777 nm triplet lines are 
arguably the most reliable since the others are either weak or partly blended. 
The main challenge when employing the O i lines relates to quantifying the non- 
LTE effects, which can be significant. The predicted non-LTE abundance cor- 
recti ons are little depend ent on the model atmosphere in question, be it ID or 
3D ( Asplund et al. 20041 ) but they do depend on the adopted coUisi o nal cross- 
sectio ns, in particular the inelastic H collisions (e.g. Asplundl 2005, Kiselman 
I993I ). Different approaches have been used in solar abundance works to date 



where some neglect the H collisions altogether based on t he available at omic 
physics data for other elements while other use the classical iDrawinI (1196811 for- 
mula, possibly with a scaling factor 5h tha t typically varies froni to 1. iHolweger 
(|200lh found logeo = 8.71 ±0.05 using the lHolweger fc Miillerl (Il974l'l model with 



granu lation corrections estimated from a 2D solar model, while lAsplund et al 



( 2OO4I ) obtained a 0.07 dex lower abundance from a 3D non-LTE study; the main 
difference, however, bet ween the two studies is the adopted H collisions (5h = 1 



:e, nowever, oer ween tne two stuaies is tne aaoptea n comsions ( jh = i 
respectively). Allende Prieto. Asplund &: Fabiani Bendichol ( 20041 ) at- 



and 0, 

tempted to calibrate the poorly known collisional efficiency through a comparison 
of the center-to-limb variation of the Oi 777 nm line strengths since the rela- 
tive importance of H collisions change as Nyi/N^ increases with smaller optical 
depths. They found that 5h = 1 was a marginally better fit than = 
while the LTE case could b e ruled out at high significance level. Recently, 
Pereira. Asplund &: Kiselman ( 20091 ) have revisited the issue equipped with bet- 



ter solar observations and tested additional 5h values but have arrived at a similar 
conclusion, as can be seen in Fig. [5l 

The Oi-based abundance we adopt here stems from Asplund et al. ( 2009d ). 
which is based on full 3D non-LTE ca lculations with S y[ = 1. The new electron 
colli sional cross-sections computed by Barklem ( 2007al ) have also been adopted 
(see Fabbian et al. 20091 for an account of the i mpact of these on ot her stars). 
The revised abundance is slightly higher than in Asplund et al.l ( 20041 ) : logeo = 
8.69 ± 0.05. Using a different 3D model, ICaffau et al. ( 2008al ) obtained logeo = 
8.73 lb 0.06 (unweighted mean) for their disk-center Ol lines, in spite of using 
5h = 1/3; with 5h = 1 their result would be 0.02 dex higher. Because they 
were restricted to ID non- LTE calculatio n s, they were unable to carry out line 
profile fitting as done by lAsplund et al.l (|2009d ). They were therefore forced 
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to rely on measuring equivalent widths, which were surprisingly large compared 
to all previous published studies. Reassuringly, when the same input d ata are 



ado pted (e.g. line strengt hs, H collisional efficiency), the Caffau et al. ( 2008al ) 



and I Asplund et al.l (|2nn9d ) analyses are in excellent agreement, which shows that 
t he exact choice of 3D model is of li t tle co nsequ ence for the permitt ed O i lines. 
Grevesse. Sauval &: van Dishoeckl (|l984l ^ and ISauval et al.l introduced 



the vibration-rotation and pure ro tation lines of OH as solar O abundance indi- 
cators. In their 3D-based analysis, Asplund et al. ( 20041 ) found that the OH lines 
implied an O abundance only slightly lower than the atomic transitions. While 
the vibration lines showed no trends with excitation potential or line strength, 
the pure rotation lines did. This is not too surprising, since these strong rotation 
lines are formed in very high atmospheric layers, which are the most challenging 
to model realistically due to for example m agnetic fields, non-LTE and non- 
equilibrium chemistry. Asplund et al.l ( 2009c ) revisited the OH lines equipped 
with a new 3D model, which yielded logeo = 8.69 and 8.69 for the vibration- 
rotation and pure rotation lines, respectively; the trend with line strength for the 
latter were still present but somewhat reduced and now in the opposite direction. 
Unfortunatel y, the OH li n es hav e not yet b e en an alysed using the alternative 
3D model of ICaffau et all (|2008al l. iMelended (|2004l l added a few vibration lines 
from the first overtone band of OH, which we however do not consider sufficiently 
reliable abundance indicators due to their extreme weakness. 

The 3D-based solar O abundance that we recommend here, logeo = 8.69 it 
0.05, is a mean of the [ Oi], O i, OH vibration-rotation and OH pure rotation 
results of lAsphmd et al.l (l2009cl ). The agreement between the different abundance 
i ndica tors is very satisfactory, which is not the case with the iHolweger &: Miiller 
(Il974l ^ model as seen in TableO As explained above, the use of a 3D model is only 
one factor in the downward revis i on of the O abundance from the value 8.83 ± 
0.06 given in Grevesse &: Sauva 3 (Il998l ^ with refined non-LTE line formation, 
better accounting of blends and improved atomic and molecular data also playing 
important roles. 



Using CO lines, 
16Q/18Q ^479_^29. 



Scott et al.l (|2006l ) determined an isotopic abundance of 



Th e sign ificantly lower, non-terrestrial ratio obtained by 



Ayres. Plymate &: Kelled (120061 ) is likely the result of their use of ID model at- 
mospheres. Several lines due to ^^C^^O are present in the solar infrared spectrum 
but due to their extreme weakness we do not believe a trustworthy ^^O/^'^O ratio 
can be inferred from them at this stage. 



3.5 Intermediate Mass Elements 



Sodium: lAsphmd etaD (|2009bl ) based their 3D solar Na abundance determina- 
tion on five weak Nal lines, which all yielded consistent results (cr = 0.03 dex). 
The predicted departures from LTE in ID are slightly negative (^i —0.04 dex) and 
agree well between different studi es ('Allend e Prieto. Asplund &: Fabiani Bendichol 



20041 . IShi. Gehren fc Zhaol l2004l . Ijakcda et_alJ |2003|) . The Na abundance we 



adopt here is slightly larger than the preliminary value presented in lAsplund. Grevesse Sz Sai 
(|2005l l. 

Magnesium: The abundance of Mg can be inferred from two ionization stages. 



each with a multitude of seemingly clean, relatively weak lines ([Lambert Sz Luck 
19781 ). Unfortunately, the (//-values for Mgi and Mgii are notoriously uncer- 



tain, as refiected in fairly large line-to- line abundance variations for both species. 
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Departures from LTE i n ID line formation by Mgi have been considered by 
Zhao. Butler &: Gehren ( 19981 ) and T. Gehren (2 008, private communi c atioii ) ; the 



abundance corrections are small and positive. lAbia k Mashonkinal (|2004l ) and 
L. Mashonkina (2008, private communication) have made corresponding calcula- 
tions for Mgii. With the predicted non-LTE effects, the mean Mgii abundance 
is brought into better agreement with the meteoritic abundance while Mgi is 
further offset. Since we do not have greater confidence in either set of transition 
probabilities and non-LTE calc ulations are available for both species, we simply 
take a straight mean of all lines ( Asplund et al. 2009bl ) . The attached uncertainty 
is disappointingly large for such an important element. 

Alu minium: Our recommended Al abundance comes from Asplund et al.l 
(|2009b ). who performed an analysis of seven good, weak All lines with relatively 
well-determined ^'/-values. The departures from LTE are predicte d to be small. 
We h ave adopted the non-LTE results from ID calculations by Gehren et al.l 
( 20041 ) ■ with the expectation that they will not be very different in 3D. 

Silicon: Silicon is a critical element since it p rovides the abso lute normaliza- 
tion for the meteoritic abundances (Sect. 14. ip . Asplund ( 2000l ) carried out an 
early 3D-based LTE analysis of Sii and Sin lines and found logesi = 7.51 ±0.04. 
With an impro v ed SD solar model and non-LT E corrections (~ —0.02 dex) taken 
from Shi et al. I (I2OO8I ). lAsplund et aD (|2009bl l found the same value in their re- 
analysis, which we adopt here. Had we instead opted to use the lHolweger &: Miillei 
(|l974l ) ] D semi-empirical model the result would have been a mere 0.02 dex 
higher. 

Phosphorus: ICaffau et all (|2007bl ) carried out a solar P determination using 
a co5bold 3D solar model and found log ep = 5 .46 (a = 0.04) from fou r P i lines. 
A slightly lower 3D-based value was obtained bv lAsphmd et al.l (l2009bl ). who also 
included one more line: logep = 5.41 it 0.03 {a = 0.03). We adopt this latter 
value on account of its smaller line-to-line scatter. No non-LTE investigation has 
been performed to date for these high-excitation P i lines. Departures from LTE 
are not expected to be significant howe yer, a prediction sup ported by the small 



200,5l'l. 



non-LTE effects for the similar S i lines (iTakeda et al.l 

Sulphur: In their 3D-based analysis, Caffau et al.l (j2007al ) obtained a dis- 



turbingly large line-to-line scatter: log es 



the non-LTE corrections of Takeda et al 



= 7.2 1 ±0.11 (a); this result includes 
(1200,511. V ery rec ently, the issue of the 



solar S abundance was revisited by Asplund et al. ( 2009bl ). who included a few 
additional lines. Their non-LTE abundance is slightly smaller and, importantly, 
shows a much re duced uncertainty: log es = 7.12 ± 0.03 (a = .03). The dif- 
ferences between Asplund et al. ( 2009bl ) and Caffau et al. (|2007al ) can be partly 
traced to the adopte d equivalent wi d ths, bu t this cannot be the whole expla- 
nation. Interestingly, Asplund et al.l ( 2009bl ) concluded that the forbidden [Si] 
1082 nm line returns too high an abundance by 0.3 dex, which may signal an 
unidentified blend, non-LTE e ff ects o r that the theoretical transition probability 
Caffau k Ludwid (|2007l l on the other hand found excellent agree- 



is m error. 



ment with the meteoritic value for this line. They adopted however an outdated 
gf -value and a much smaller equivalent width, which we are unable to reconcile 
with the observ ed solar spectrum. 

Potassium: Asplund et al.l ( 2009bl ) based their 3D analysis on six Ki lines, 
including the strong resonance line at 769.9 nm. As in previous studies, the line- 
to-line scatter is disappointingly large for these lin es, even when taki ng into ac- 
count the theoretical non- LTE effects computed bv lZhang et"all (|2006l l. In terms 
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of abundances, the non-LTE corrections are substantial for the two strongest hues 
(> -0.2 dex). 

Calcium: The solar Ca abundance can be determined from both Ca i and Ca ii 
( Lambert Sz Luck 19781 ). although there are problems with both species: most of 
the Cal lines are rather strong, whilst the uncertainties in the transition proba- 
bilities for the Call li nes are still uii s atisfact ory large (Opacity Project data as- 
suming LS-coupling). lAsplund et al. found good agreement between the 
two ionization stages and obtained a final Ca abundance of logeca = 6.34 it 0.04 
(cj = 0.04) from all C al and Call lines, includin g ID non-LTE abundance cor- 
rections computed by Mashonkina. Korn Sz Przybillal (|2007i ). We note that the 
forbidden [Can] line at 732.3 nm implies an abundance which is fully consistent 
with the other transitions. 



3.6 Iron-peak Elements 



Scandium: With its relatively low ionization potential, Sc is prone to over- 
ionization via non-LTE effects ( Zhang. Gehren Zhao 20081 ). The impact on 
abundances from Sc i lines is dramatic (-1-0.15 dex relative to LTE on average 
for the lines considered by Scott et al. 2009bl ) . but predictably mild for Sen 
(—0.01 dex). Good lines, oscillator strengths and hyperfine constants are available 
for both ionization sta ges. With five weak lines of S c i, twelve of Sen an d the 
non-LTE corrections of lZhang. Gehren fc Zhaol (|2008l ^. IScott et al.l (l2009bl ^ have 
produced 3D-based abundances showing very good agreement between the two 
ionization stages. In Table [1] we give the mean of the abundances from all lines of 
both species, where Sc ii ultimately dominates because of the greater number of 
lines. The -1-0.10 dex difference between the new Sc abundance and the meteoritic 
one is intriguing. 

Titanium: Very many good lines of Til are present in the Sun, virtually 
all with excellent atomic data. Whilst fewer lines of Tin are available and the 
atomic data are not quite as good as for Til, the situation is still relatively 
acceptable. The one glaring absence is a study of non-LTE effects for this element. 
As is the case for its neighbours Sc and V, the low ionization potential of Ti 
and the resulting minority status of Til are expected to significantly skew the 
abundances it retu rns under t h e assu mption of LTE. Using the improved 3D 
model atmosphere, Scott et al.l ( 2009bl ) found stark disagreement between LTE 
abundances from Til and Tin (Til: logexi = 4.85 ± 0.06, Tin: logexi = 4.99 ± 
0.04), though with very low scatter from each species individually (cr = 0.03 for 
both stages). In the absence of any non-LTE study, we favour Tin in Tabled] by 
taking the mean of the two results weighted according to their total uncertainties; 
we only retain any weighting for the Til result because of the sheer number of 
excellent lines and oscillator strengths it includes. 

Vanadium: Neutral vanadium has scores of excellent lines in the visible solar 
spectrum, whilst very few are available from V ii, all extremely perturbed. On the 
other hand, Vl is expected to show similarly strong non-LTE effects as verified 
for Sci and suspected for Til, owing to the elements' similar potential for over- 
ionization. The impact of no n-LTE line f ormati on on V has not yet been the 
subject of a dedicated study. Scott et al. ( 2009*9 ) obtain a 3D LTE abundance 
of logev = 3.83 from Vi lines. Scott et al.l 's analysis of Vii lines yields an 
abundance quite consistent with this corrected value, but the poor quality of the 
lines employed makes the result inconclusive. In the absence of any good results 
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for Vii, we adopt the Vi value but adjust it by the mean non-LTE effect seen in 
Sci and Cri (+0.10 dex). 

Chromium: The solar spectrum exhibits a large number of clean, weak lines 
of both neutral and once-ionised chromium. Unfortunately, no reliable oscilla- 
tor strengths are available for Cm. Neither is any estimate of the impact of 
non-LTE effects on Cr available; given its relatively low ionization potential, Cr 
should consist mostly of Cr ii in the Sun, possibly leading to non-LTE effects for 
Cri. Recent non-LTE calculations by M. Bergemann (private communication) 
confirm departures from LT E at the level of -1-0.04 dex. Using their own measured 
oscillator strengths for Cri. lSobeck. Lawler &: SnedenI (1200711 recently revised the 



(|2009bl l 



US- 



solar chromium abundance. This result was updated by IScott et al. 
ing a 3D model. Given the much greater number of Cri lines in the sample, and 
the poor quality of the Cm oscillator strengths, we allow Cri to dominate the 
value in Table [J by taking the mean of abundances from all Cr i and Cr ii lines 
together. The only reason we allow abundances from Cr ii lines to enter into the 
mean is our suspicion of an as yet unquantified non-LTE effect for Cri. 

Manganese: Following years of confusion now traceable mostly to poor os- 
cillator strengths, the solar manganese a bundance was recently rederived by 
Blackwell- Whitehead &: Bergemann (|2007l l. using their own improved Mnl q f- 



values , extensive hyperfine splitting data and the non-LT E calculations of lBergemann fc Gehren 
(|2007l ). This result was updated bv lScott et all (|2009bl l using a 3D solar model, 
bringing the photospheric Mn abundance even closer to the meteoritic value. Mn i 
shows moderate non-LTE effects (about +0.06 dex), consistent with its ionization 
potential being intermediate between the lower values seen in the lighter mem- 
bers of the iron group Sc, Ti, V and Cr and the higher values exhibited by the 
heavier elements Fe, Co and Ni. 

Iron: In many circumstances iron serves as a proxy for the overall metal 
content and as a reference against which other elemental abundances are mea- 
sured. An intense debate long raged as to whe ther the solar Fe abundance in- 
ferred from Fe I lines should be logepp ~ 7.7 (e.g. iBlackwell. Lynas-Gray &: Smith 
19951 ) or logepe ~ 7.5 (e.g. iHolweger. Kock fc Bardlll99,^ l. with the choice of 



equivalent widths, (//-values, microturbulence parameters and pressure-damping 
constants all play ing a role. The adve r it of improved t ransition probabilities 
for Fell lines (e.g. Hannaford et al. 19921 . Schnabel. Kock & : Holweger 1999) and 
a proper treatment of pressure broa dening (e.g. Anstee. Q'Mara &: Rossi 1997 . 
Barklem. Piskunov fc O'Marall2000"al ) strengthened the case for a low Fe abun- 



dance. An unsatisfactory combination of large abundance scatter, the presence of 
a free parameter (microturbulence) , and a t rend with excitati on potential still re- 
mained with the iHolweger fc Miillei] (|l974l l model, however (iGrevesse fc Sauvall 
19991'). Using a 3D h ydrodynamical solar model and LTE line profile fitting. 



Asplund et al.l (l2000d ) succeeded in achieving consistent results for Fe i and Fe ii 

lines without invoking any microturbulence: log epe = 7.45 ± 0.05. 

The solar Fe abundance we advocate here comes from Scott et al. (l2009bl l. 
The Fei-based abundance of l og epe ^ 7.52 ± 0.05 from weak lines is slightly 
higher than in Asplund et al. ( 2000a ) on account of a slightly shallower mean 
temperature structure in the new 3D model and allowance of departures from 
LTE. T he non-LTE abundance corrections have been estimated to be +0.03 dex 
for Fe I (|Collet. Asplund Theveninll2005l : iKorn. Shi GehrenI 120031') . although 
this depends on the still poorly-known inelastic H collisions (see Asplund 20051 . 
and references therein). In fact, a more reliable abundance now comes from 
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Fell li nes, thanks to the improved transi tion probabihties of Melendez Sz Barbuvl 
and the broadening treatment of lBarklem &: Aspelund- Johansson! (|2005l ). 
Together these imply l ogepe ^ 7.50 ifc O.OMa = 0.04) with the same 3D model. 
We note that with the Holweger Sz Miilleil ( 19741 ) model, there is a trend in Fei 



abundances with excitation potential and a discrepancy of 0.14 dex between Fei 
and Fell. Both have disappeared with the 3D solar model (Fig. [6]). Reassuringly, 
the situation regarding Fe has greatly improved over the past decade. 

Cobalt: The fairly high ionization potential of Co and the availability of 
many good Col lines, q.f-valu es and hype r fine co nstants should make its solar 
analysfs quite straightfo^d. IScottetal.1 M ) found a 3D LTE abundance 
of log eco = 4.88, in agreement with the m eteoritic value. Ta king into account 
the surprisingly large non-LTE corrections of lBergemann ( 20081 ) changes the final 
abundance to log eco = 4.99 ± 0.07 (a = 0.05). 

Nickel: The solar Ni abundance was estimated bv lScott et aD (|2009al ) from 17 
weak Ni i lines, showing a very sm all abundance scatt er. The value given in Table 
[T]is a slight revision of this result ( Scott et al. 2009bl ) based on the improved 3D 
model atmosphere that has been employed throughout this review. 



3.7 Neutron Capture Elements 

Quite some work on the solar abundances of the heavy elements has been car- 
ried out over the past decade, largely driven by improvements in atomic tran- 
sition probabilities. In this context, we would like to particularly highlight the 
very careful labora tory measurements performed by the Wisconsin group (e.g. 
Lawler et al]l2009l ). In a beautiful series of papers, they have systematically 
and methodically measured radiative lifetimes, branching fractions, oscillator 
stre ngths, isotopic shifts and hyperfin e con stants for all the ra re E arth elements 



La (Lawler. Bonvallet &: Sn eden"200ll). Ce ( Lawler et a 



20091) Nd (iDen Hartog et al.. .20031) . Sm (iLawler et al 
2001bl ). G d (iDen Hartog et al.ll20oi) Tb (iLawler et al.l 



20091 ). Pr fSneden et al 



20061). Eu (Lawler et al. 



2001d). Dv (jSneden et al 



20091). Ho (ILawler Sneden fc Cowanll2004l ). Er (ILawler et ah l boosl ). Tm, Yb and 
Lu (ISneden et al.1 120091 )7 the group has also carried out a similar study for Hf 



(ILawler et al.l l2007l)~ In t h eir so lar abundance analyses, they have employed 



the ID Holweger fc Miilleil (|l974l ) model atmosphere and assumed LTE for the 
spectrum synthesis. Although all the utilized transitions are from the majority 
ionization stage (once ionized) for the rare Earths, one still expects a minor tem- 
perature and thus model-atmosphere dependence, since the lines originate from 
atomic levels with low excitation potentials (xcxc ^ 1.5 eV). The recommended 
solar abundances provided in Table [Tj for these elements have been modified 
from the original sources by calculating the abundance difference betwe en a 3D 
hydro dynamical model atmosphere and the Holweger &: Miillei ( 19741 ) model, 
each computed using the same line formation code with consistent input data 
( Grevesse et al.r2009l ). For the majority of these lines, such (3D-HM) abundance 
corrections a mount to ~ —0.04 dex. I n add ition, the ID non-LTE effects on Euii 
predicted by Mashonkina &: Gehren ( 2000l ) have been taken into account in the 

3D result; non-LTE studies are not available for an y other rare Earth elements. 

A s imil a.r procedure has been impl emen ted for Cd (lYoussef. Doenszelmann &: Grevesse 
1990|) . W (iHolweger Werneil ll983). Os (jOuinet et al.ll2006l ). Ir (|Youssef k Khahl 
198a ) and Au (|Youssejll986l ). all of which have only heavily blended transitions 
available in the solar spectrum, necessitating careful and time-consuming spec- 
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trum synthe s es. D ue to blends, we consider the two Asi hnes employed by 
Gopka et alJ (120011 ) too unreliable to base a meaningful abundance determina- 
tion on and consequently do not provide a solar photospheric As abundance 
here. Similarly, we are unable to defend an estimate of the Sb abundance from 
th e very weak and hea vily perturbed, purported Sbi line at 323.2 nm reported 
bv lRoss fc AUctI (Il97fil ). 

For all other neutron capture elements, the adopted abundances come from 
a complete 3D-based reanalysis using the best atomic data and newly measured 



equiva lent widths or spectrum synthesis for the most reliable lines. iGrevesse et al 



( 20091 ) present the solar analyses of Cu, Zn, Ga, Ge, Rb, Sr, Y, Zr, Nb, Mo, Ru, 



Rh, Pd, Ag, Cd, Sn, Ba, Pb anc 


Th, which we adopt throughout; the case of 


Ru is also discussed in Fivet et al 


(20091). The available ID non-LTE abundance 


corrections for Znl (Takeda et al. 


I2OO5I). Srii (Mashonkina & Gehrenll200ll) and 


Ball (Mashonkina &: Gehren 200d) have been included in these estimates. For 



comparison purposes we have also repeated the line calculations with an assort- 
ment of ID model atmosphere s. We have found very good agreement (< 0.02 dex) 
with recent studies using the iHolweger k Miilled (|l974l ) model when allowance 
is made for differences in ad opte d transition probabilities a nd e quivalent width s 
for Ge (iBiemont et al.lll99 9l). Sr (iBarklem Q'Maral[200ol ). Pd (|Xu et al.l[20oi ) 
and Pb (|Biemont et al.l l2000). 

A handful of elements hav e prev iously been exposed to a 3D-based solar abun- 
dance analysis. Ljung et al.l ( 20061 ) made use of new tran s ition p robabilities for 
Zrii and employed the 3 D solar model o f lAsplund et al. ( 2000bl ). Not surpris- 
ingly, t he agreement with iGreyesse et ah ( 20091 ) is very good. Eu has been the 
focus of lMucciarelh et al.1 (|200 8'). who employed a 3D co5bold hydro dynamical 
model and found a mean Eu abundance iii perfect agreement with ours. The 3D 
LTE abundance for Hf fro m iGaffau et al.1 (l2008bl ) is only 0.02 dex higher than 
the corresponding value by Grevesse et al.l ( 20091 ) . who corrected the spectrum 
synthesis results of Lawler et al. ( 200?! ) for 3D effects rather than at tempting to 



measur e the equivalent widths of these weak and blended Hf 11 lines. iCaffau et al. 



(|2008bl ) also performed a detailed investigation using 3D line profile fitting of the 
Thii 401.9 nm line, which is distorted by two Coi and Vi lines (not Vii as stated 
in their paper), as well as being in the red wing of a relatively strong feature 
from Fei and Nil. Interestingly, they argued that it is important to take into 
account the l ine asymmetry of the Fel+Nii blend introduced by the convective 
motions (e.g. Asplund et al.l 2000bl ). which otherwise woul d lead to an o yeresti - 
mation of the Th abundance by ~ 0.1 dex. Unfortunately Caffau et al. ( 2008bl ) 
did not specify what abundances they used for the blend ing lines, so it i s diffic ult 
to exactly reproduce t heir result; th e value obtained by Grevesse et al. ( 20091 ) is 
consistent with that of lCaffau et al.l . though smaller by 0.06 dex. 



3.8 Abundances from Sunspots: F, CI, In, Tl 

The photospheric abundances of a few elements are not possible to measure us- 
ing spectroscopy of the quiet Sun. The solar fluorine and chlorine abundances 
have in stead been iiiferred from the i nfrared spectrum o f sunspots using lines 
of HF (iHall k Noveslll969l ) and HCl (jHah k Novedll972l ). Given the improve- 
ments in molecular data, observational quality and sunspot modelling in the 
intervening years since the original workds, it seems like a critical re-evaluation 
of these studies is warranted but not attempted here. The same holds for Tl 
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(|Lambert. Mallia k Sm"ithlll972l l. We note that in the meantime a more rehable 
estimate of the solar CI may be pro vided by the abundance of nea rby Hii re- 
gions, whic h is lo g eci = 5.32 ib 0.07 ( Garcia- Rojas &: EstebanI 2007 ). Recently, 
Vita.eta,l.l have revisited the solarin abundance. In contrast to previous 

studies based on the quiet Sun, they find an abundance from a sunspot spectrum 
in agreement with the meteoritic value. 



3.9 Indirect Determinations: Noble Gases 

Due to their high excitation potentials, there are no photospheric lines of the noble 
gases helium, neon, argon, krypton and xenon available in the solar spectrum, 
which forces one to estimate the solar abundances of these elements through 
indirect methods. 

Helium: A highly accurate He abundance determination is possible through 
helioseismology. The technique utilizes the change in the adiabatic index Fi in the 
He II ionization zone that takes place at a solar radius of r/i?© ~ 0.98. The result 
is quite insensitive to the reference solar mode l employed for the in version, but 
does depend on the adopted equation-of-state. Basu &: Antial ( 20041 ) found a He 
mass fraction in the convection zone of Kj = 0.2485 ± 0.0034, which corresponds 
to a He abundance by number of logenc = 10.93 ifc 0.01. The uncertainty es- 
senti ally reflects t he difference between the MHD ( Mihalas. Dappen fc Hummed 
19881 ) and OPAL ( Rogers &: Navfonovll2002l ) equation-of-state, but there is also 
a minor depe ndence on the assumed metallicity of t he reference solar model for 
the inversion ( Trampedach. Dappen &: Baturinll2006l ) . A similar value is obtained 
when calibrating solar interior models to achieve the correct m ass, luminosity and 



radius at the present solar age, if the chemical composition of lGrevesse &: Sauval 
(119981') is adopted and diffu s ion is considered. With the solar abundances of 



Asplund. Grevesse &: Sauval ( 20051 ) on the other hand, the calibrated proto-solar 
He abundance would imply a surface value lower than the helioseismic measure- 
ment by 0.02 dex. Such models are also in conflict with other helioseismological 
evidence such as the depth of the convection zone and the sound speed variation 
with depth, as discussed in Sect. 14.31 

Neon: The Ne abundance can be inferred from X-ray and UV spectroscopy of 
the solar corona and solar flares as well as directly from the solar wind. The inter- 
pretation is, however, complicated by the so-called first ionization potential (FTP) 
effect: elements with ionization potential Xion < 10 eV are enhanced in the upper 
solar atmosphere and solar wind compared with their photos pheric values, whilst 
the h igh ionization elements are not or only slightly affected ( Feldman &: WidingI 
2003). Rather than measuring the Ne abundance directly, a common approach 
is to determine the abundance relative to a reference element and assume that 
the ratio is the same in the photosphere. Here O is the reference element of 
choice, although we note that there is still a significant difference in ionization 
potential between O (13.6 eV) and Ne (21.6 eV), which could potentially intro- 
duce a systematic error given our poor understanding of the physical reasons for 
the FIP effect. The degree of chemical sepa ration varies sig nificantly, being more 

found a Ne/0 ratio 



Yound (|2005 



severe in regions of higher solar activity, 
of 0.175 ± 0.031 for the quiet Sun, which is the value we adopt here to convert 
our photospheric O abundan ce to a Ne con tent. This result is particularly ap- 
pealing for our purposes since I Young (|2005al l has demonstrated that the average 
quiet Sun does not show any FIP effect, a conclusion which is further corrob- 
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orated when using the solar abundance s given herein. This ratio is consistent 
with the results of ISchmelz et al" I (I2OO5II for active regions (0.175 it 0.074) and 
for the solar wind (0.14 ± 0.03) bv iBochslerl (|2007l l. The corresponding ratio in 
solar ene rgetic particles is very similar and remarkably constant: 0. 1 52 ± 0.006 
( Reamed il999). which was adopted by Asplund. Grevesse Sauva 1 (|200,5l V In 
spite of its very small uncertainty, here we weight this result less than from the 
quiet Sun due to the potential modification through the FIP effect. Together 
with our recommended O abundance, the photospheric Ne abundance becomes 
logeNe = 7.93 ±0.09. 

Our value is smaller than t he sol ar flare determination (logeNe = 8.11±0.12) by 
Landi. Feldman k Dosch^ (120071 1. but we consider the quiet Sun a more reliable 



indicator than solar flares due to possible fraction ation processes in the latter. 
While the Ne abundance (logeNe = 7.96 ± 0.13) of IBochslerl (|2007l l agrees with 
ours, there are significant additional uncertainties in this result on account of 
using He as a reference element, as it is known to be highly variable in the solar 
w ind. 

Drake &: Testa! (|2005l l have argued on the basis of X-ray spectroscopy of nearby 



stars that the Ne/0 ratio is significantly higher (0.41) than the solar ratios dis- 
cussed above. Their sample is, however, heavily biased towards active stars, 
which experience an inverse FIP effect: elements with Xion > 1 eV are enha nced 
in th e corona compared with low i onizat ion potential species (|Giidelll2004l l. In- 
deed, iRobrade. Schmitt &: Favata ( 20081 ) show that there is a systematic trend 
in Ne/0 with activity level such that other inactive stars have similar low ratios 
to the Su n. One would t herefo re suspect that the high coronal Ne/0 ratio in- 
ferred by iDrake &: Testal (120051 ) is a consequence of the difference in ionization 
potential, making Ne more susceptible to fractionation than O in high activity 
environments. As will be further discussed in Sect. 14.21 & more reliable stellar 
gauge of the solar Ne abundance is provided by OB stars in the solar neigh- 
borhood. Both the Ne/0 ratio (0.20 ± 0.03) and the absolute Ne abundance 
(logeNe = 8.08 ±0.03) estimated in this way are in very good agreement with the 
solar values we adopt here when the effects of solar dif fusion and Galactic chemi- 
cal en richment over the past 4.5 Gyr are accounted for (jPrzybilla. Nieva &: Butler 



20081). 



Argon: As discussed in detail by Lodders ( 20081 ). there are a wide range of 
methods for inferring the solar Ar abundance: solar wind measurements, solar 
flares and energetic particles, nuclear statistical equilibrium with the abundances 
of the nearby a-elements ^^S and ^'^Ca, and comparison with other solar system 
objects (Jupiter) and the solar neigh borhood ( B sta rs, planetary nebulae, Hii 
regions). We follow the procedure of Lodderd ( 20081 ) but employ our preferred 
values for the photospheric abundances of O (for the solar wind, flares and ener- 
getic particles), Kr and Xe (for analyses of Jupiter), and the lower He abundance 
appropriate for the upper atmosphere instead of the photospheric value (for solar 
wind analyses). The quoted errors for all methods in the literature are only the 
dispersions of the measurements, so likely underestimate of the true uncertain- 
ties. We therefore take a straight unweighted mean of all techniques, wh ich leads 
to log eAr = 6.40 ± 0.13. Our value is 0.1 dex lower than th e resul t of Lodderd 
(I2OO8I ) but 0.2 dex higher than lAsplund. Grevesse k Sau val' (200^ , which was 
based only on the Ar/0 ratio in solar energetic particles (,Reameaill999l ). 

Krypton: The Kr abundance has been estimated from interpolation of the the- 
oretical s-process production rates, since the cross-sections for neutron capture 
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of nearby pure s-isotopes have been accurately measured experimentally (Palme 
& Beer 1993). Our value accounts for a slight revision of the photospheric Si 
abundance, which serves as a normalization of the s-process predictions. This Kr 
estimate agrees well with the value inferred from the solar wind Kr/Q ratio to- 
gether with our ph otospheric O abundance (iGeiss. Gloeckler &: von Steigeiill994l ; 
Heber et al.ll2nn8bl ). 

Xenon: The procedure to estimate the Xe content is the same as for Kr 



(lKappeleiil2005l ) . For Xe, however, this s-process estimated value i s 0.4 dex lower 

than deduced from the solar wind as measured from luna r regoliths (jGeiss. Gloeckler &: von i 
19941 ) and the Genesis space probe (|Heber et al.ll2008bl ). This may suggest that 
Xe has experienced fractionation through the FTP effect in spite of its relatively 
high ionization energy (xion = 12.1 eV). 



3.10 Solar Isotopic Abundances 

As a convenience to the reader, in Table [3] we give our best estimates of the proto- 
solar isotopic fractions for each element. In most cases the terrestrial isotopic 
abundance ratios as recommended by International Union of Pure and Applied 
Chemistry ( Rosman fc Tavlor 19981 ) have been adopted as representative of the 
solar system values; direct solar isotopic information can only be g leaned from 



CO v ibration-rotation lines in the infrared spectrum of the Sun (e.g. IScott et al. 



20061 ) . although with comparatively large error bars by terrestrial standards. It 
should be noted that except for C and O, the isotopic fractions given here are the 
more conservative "repre sentative isoto pic composition " from the compilation of 
Rosman &: Tavlor ( 19981 ). in contrast to Lodders ( 20031 ) who opted for the "best 
measurement from a single terrestrial source" from the same reference. For most 
volatile elements more reliable isotopic fractions are available from alternative 
solar system sources due to strong depletion of terrestrial samples, as discussed 
below. 

The proto-solar deuterium abundance can be estimated from observations of 
Jupiter. Using ISO spectra of the Jovian atmosphere and t aking into account an 
expected 5 — 10% enrichment during the planet's formation. iLellouch et al~ ( 200ll ) 
obtained a proto-solar abundance of D/H = (2.1 ± 0.4) ■ 10~^. An alternative 
approach makes use of the measured enrichment of the ^He/^He ratio in the 
solar wind in comparison w ith the proto-solar r atio ^He/^He = (1.66 it 0.05) • 
10~^ inferred from Jupiter ( Mahaffy et al.^^1998 ), since the ad ditional ^He can 
be attributed to D-burning in the Sun (Gloeck ler fc Geisd l200nl ) . Together with 
the solar wind ratio ^He/'^He = (4.53 it 0.03) • 10~^ from Genesis measurements 
(jHeber et all l2008al ). the Jovian ratio implies D/H = (1.96 ± 0.3) • 10"^; most 
of the uncertainty comes from the poorly-known He fractionation in the solar 
wind. From the two methods, we finally estimate a proto-solar abundance of 
D/H = (2.0 lb 0.2) ■ 10 ~^. This is consiste nt with the present-day local interstellar 
medium abundance ( Linskv et al.l 20061 ). which implies only a small degree of 
astration over the past 4.5 Gyr. 

The ^^C/i^c ratio does not appear to vary significantly between different solar 
system sources, so we simply adopt the accura te terrestrial referen ce value of Vi- 
enna Peeddee belemnite: 12C/13C = 89.4±0.2 (ICoplen et al.ll2002l'l We no te that 
the "representative isotopic composition" from Rosman &: Tavlor ( 19981 ) corre- 
sponds to the significantly more uncertain value ^-^C/^^C = 92 it 7. The ratio 
inferred from the solar photosphere using CO lines (86.8 it 3.8) is consistent with 
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the terrestrial value ( Scott et al. 20061 ) . 

The ^^N/^^N r atio varies tremen dously between various solar system objects, 



from 100 to 450 (jOwen et al.ll200ll ). Here we adopt the Jovian value 435 it 57 as 



representative of the proto-solar ratio, in agr eement with va lues in the interstellar 
medium and carbonaceous chondrites ((Meibom et al. 20071 ). 

Oxygen isotopic abundances display a variability of up to a few percent within 
the solar system. The best val ues to adopt for the proto-solar oxygen iso- 



topic ratios are still debated (e.g. iHashizume &: ChaussidonI l2005l . [Ireland et al 
20061). The first resu lts from the Genesis mission should appear very soon 



( McKeegan et al. 20081 ). which will hopefully shed some light on the issue. Until 
such further clarification, we adopt the terrestrial values represented in Vienna 
Standard Mean Oc ean Water: ^^O/^^O = 498.7 ± 0.1 and ^^O/^^O = 2632 ± 7 
(ICoplen et al.ll2002ll. The solar photospheric ratio ^^O/^^^O = 479 ± 29 obtained 
by Scott et al.l ~( 20061 ) is consistent with the terrestrial value, but with larger 
uncertainties ( although in fairnes s the t errestrial "representative isotopic com- 
position" from iRosman &: Taylor (Il998l) is equally uncertain: 487 it 35). Sur- 
prisingly, Ayres. Plvmate &: Keller ( 20061 ) found strongly sub-terrestrial ratios 
(e.g. ^^O/^^O = 440 lb 20) from their analysis of the solar CO lines, probably 
due to an inadequate consideration of the differing ter nperature sensitiv ities of 
the isotopomer transitions with ID model atmospheres ( Scott et al. 20061 ). Lines 
from ^^C^^O are also visible in the solar infrared spectrum, but due to their 

extreme weakness their reliability as abundance indicators is ques tionable . 

For Ne, Ar, Kr and Xe we adopt the recommended values of IWieleil (|2002l ) 
derived from solar wind and lunar sample analyses. Very similar ratios have re- 
cently been pres ented from prelimin ary studies of solar wind samples from the 



Genesis mis sion (Heber et al. I2p08b ). For Mo, Dy, Er, Yb and Lu we have ex- 



changed thelRosman Sz Tavlorl (119981) v alues with the more recent measurements 



Lodders. Palme &: Gaill ( 20091 ). Similarly the proto-solar values given 



given m 

here for the radioactive nuclides have been taken from the same source. 



3.11 Conversion From Photospheric to Bulk Abundances 

The elemental abundances given in Table [1] are those currently present in the 
solar photosphere and by extension the whole convection zone. Due to the 
combined effects of thermal diffusion, gravitational settling and radiative ac- 
celeration (often referred to collectively as diffusion) over the past 4.56 Gyr, 
the values in Table [T] differ slightly from the proto-solar chemical composition. 
In comparison with present-day photospheric values, the bulk composition of 
the Sun is higher by 0.05 dex for helium and 0.04 dex for the heavier elements 
(jTurcotte &: Wimmer-Schweingrub e3[200i). The uncertainty in the overall diffu- 
sion correction is judged to be about 0.01 dex. The element-to-element variation 
in the expected diffusion is very small (< 0.01 dex) and therefore currently not 
observationally testable by a cq r nparis on of photospheric and meteoritic abun- 
dances. We note that Lodders ( 20031') adopted a somewha t large r correction, 
0.074 dex, than that of Turcotte &: Wimmer-Schweingrubei] (I2OO2I ). She based 
her estimate on the results of Boothrovd &: Sackmann ( 20031 ). who computed 
their solar models using a less detailed diffusion model. 
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3.12 Solar metallicity 

With the solar chemical composition we recommend here, the mass fractions of H, 
He and metals in the present-day photosphere become X = 0.7381, Y = 0.2485 
and Z = 0.0134. The respective values for the bulk composition are 0.7154, 0.2703 
and 0.0142. In ot her words, the solar metallicity is no longer the canonical 2% 
recommended by Anders &: Grevessel ( 19891 ) but rather a substantially smaller 



1.4%. Indeed, the intervening years have seen a steady decrease in Z/X from 
0.027 to our preferred value of 0.018. This is illustrated in Tabled! which collects 
the resulting mass fractions from a number of widely used compilations of the 
solar chemical composition over the past two decades. Our metallicity reverses 
this declining trend of Z wi t h tirn e by being slightly larger than advocated by 
Asplund. Grevesse &: Sauva 3 Hooi), which is mainly the result of the somewhat 



higher C, O, Ne and Fe abundances. 

For com pleteness, we note that th e He abundances recommended by Anders Sz Grevessel 



(Il989l ^ and lGrevesse fc Noels! (|l993l l were based on calibration of solar models and 



hence corresponded to bulk composition; here we have substituted those with the 
helioseismic convection zone He content (Sect. 13. 9p when computing the present- 
day photospheric values . Finally , frorn Table H! it is clear that the surface He 
abundance employed by Lodderd ( 20031 ) is too low on account of adopting a too 



efficient elemental diffusion, as discussed in Sect. 13.111 

4 OTHER METHODS FOR INFERRING SOLAR ABUNDANCES 
4.1 Meteorites 

Meteoritic data on the elemental abundances constitute an excellent comparator 
for the results from the solar photosphere presented above. Not all meteorites are 
suitable for the purpose though; only a very small subsample are believed to ac- 
curately reflect the chemical composition of the proto-solar nebula, since the ma- 
jority have undergone varying degrees of fractionation and alteration. The most 
primitive, undifferentiated meteorites in this respect are the so-called CI (or CI) 
carbonaceous chondrites, of which there a re only five kn own: the Alais, Ivuna, 
Orgueil, Revelstoke and Tonk meteorites (jLoddersI [2OO3I 1. Of these, Orgueil is 



the by far largest and has therefore been exposed to the greatest number of 
analyses. It should be noted that while CI chondrites are believed not to have 
been noticeably chemically fractionated, they have still been severely modified 
mineralogically by aqueous alteration and thermal metamorphism. The main ad- 
vantage of meteorites as a cosmic abundance standard is the extraordinarily high 
precision obtained, even for isotopic abundances, through mass spectroscopy. It 
comes at a cost, however, since the volatile elements H, C, N, O and the noble 
gases - which also happen to be the most abundant elements - have been severely 
depleted. As a consequence, the abundances in meteorites are pegged to the Si 
content, rather than H as is the case in solar and stellar spectroscopy. 

The mean abundances of CI chondrites prese nted in Table U have b e en ap - 
propriated from the data carefully compiled by iLodders. Palme Gail! (Hooi). 



To place the meteoritic abundances (defined as A^si = 10^) on the same abso- 
lute scale as for the photosphere (defined as log en = 12.00), we have renormal- 
ized their data such that the meteoritic and photospheric Si abundances agree. 
The relation between the two scales is thus specified as log ex = 1.51 + log A^x- 
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In principle the conyersio n factor could employ additional elements besides Si; 
Anders k Grevessd (Il989ll for example used Nia., Mg , Si, Ca, V, Cr, Co, Ni, Y, 
Zr, Nb and Mo, while iLodders. Palme Sz Gaill (|2009l ) included as many as 39 el- 
ements. In practice, however, this does not change the result noticeably but it 
introduces some degree of arbitrariness in the selection of elements to include in 
this exercise. Fig. [7] shows the difference between the thus obtained photospheric 
and CI chondrite abundances, highlighting the excellent correspondence between 
the two methods. Naturally all volatile elements are depleted in the meteorites, 
whereas Li is depleted in the solar photosphere. For the other 45 elements with a 
purported photospheric uncertainty less than 25%, the mean difference between 
photospheric and meteoritic abundances is 0.00 it 0.05 dex. The photospheric 
and meteoritic abundances differ by more than the total combined uncertainty 
for only 10 out of the 57 non- volatile elements (expected ~ 18/57), which sug- 
gests that our estimated errors are in fact slightly too conservative. It should be 
emphasized that this excellent agreement cannot be interpreted as an absence of 
any significant diffusion in the Sun (Sect. 13. lip , as the meteoritic scale is set by 
Si; all elements heavier than He are fractionated by ~ 0.04 dex relative to H but 
hardly at all relative to each other. 

In spite of the excellent agreement overall illustrated in Fig. [TJ there are a 
few noticeable exceptions. Of the ten non-depleted elements where the solar and 
meteoritic abundances differ by more than 0.1 dex (25%), this can in all cases 
be blamed on uncertain photospheric analyses. Indeed, in many cases the two 
methods still agree marginally within the combined errors (CI, Au, Tl). We note 
a significant improvement for W compared wi th previous compilations of the solar 
chemical composition (A log e = +0.49 dex in lAsplund. Grevesse &: Sauva lHooi) 
while here the difference amount to +0.20 dex); we suspect that unidentified 
blends and/or erroneous continuum placement may be at fault for the persisting 
discrepancy, a s well as for Rb a nd Hf. For Co we reckon that the non-LTE 
corrections by Bergemann ( 20081 ) are somewhat overestimated. The remaining 
three elements Rh, Ag and Pb have had their photospheric values determined 
from transitions originating either in the ground state or from low-excitation 
levels of the neutral species (also true for Pd, W and Au), w hich are prone to 
non-LTE effects in their excitation and ionization populations ( Asplund 20051 ). 



4.2 Solar Neighborhood 

Independent estimates of the chemical composition of the Sun are provided 
by solar-type stars, OB stars, Hii regions like the Orion nebula and the in- 
terstellar medium in the solar neighborhood. In order to make a meaning- 
ful comparison one must first account for the effects of diffusion in the Sun 
and Galactic chemical enrichment. The form er amounts to 0.04 dex for met- 
als ( Turcotte Wimmer-Schweingruber 20021 . cf. Sect. 13. lip while the latter is 
typically predicted to be 0.05 — 0.15 dex dep ending on the element in question 
(jChiappini. Romano Matteuccill2003l : IPran tzos 200^). 

Whether the Sun is somehow special in its chemical composition relative to 
othe r FG dwarfs in th e Galactic thin disk has attracted attention over the years 
(e.g. Gustafsson 2008 . Robles et al. 20081 . and references therein). It should 



be borne in mind that most large-scale a. b undance analyses of Galactic thin 
disk s tars to date (e.g. Bensbv et al. 20051 . Edvardsson et al. I ll993l . [rA rmanni 
20081 . iRamfrez. Allende Prieto fc LambertI l2007l . iReddv et al.l l2003l l have been 
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performed differentially relative to the Sun to remove the uncertainties from 
^/-values and to minimize systematic errors due to inadequacies in the stellar 
modelling. These studies therefore reveal little about the absolute abundances of 
the stars but do imply that the Sun is a rather ordinary thin disk G dwarf for its 
age. In his volume-complete sample of stars within 25 pc, Fuhrmann ( 20081 ) found 

a mean spectroscopic metallicity of [Fe/H]= —0.02 ± 0.18. Similarly, among the 

stars w ithin 40 kpc having an isochrone age of 4-6 Gyr, iHolmberg. Nordstrom &: Andersen 
(I2OO8I ) obtaine d a mean metall i city o f [Fe/H]= -0.00 ± 0.10. In their study of 
11 solar twins, Melendez et al. ( 20091 ) found that the solar X/Fe ratios did not 
differ by more than 0.06 dex from the mean of the solar twin sample for any of 
the 23 investigated elements. 

In the past, the Sun appeared to be metal-rich compared with, for example, 
young OB stars in the solar neighborhood, despite the expected chemical enrich- 
ment over the past 4.5 Gyr. In the meantime, the solar abundances have decreased 
and the hot star abundances have increased. The latter have also become more 



accurate, as reflected in greatly diminishe d uncertainties (e.g. iDaflon &: Cunha 



20041 . iMorell |2008| . iNieva Przvbillal I2OO6I) . Arguably th e most accurate result 



to date is the analysis of Przvbilla. Nieva &: Butler ( 20081 ). who have carried out 
a non-LTE study of six nearby B-type main sequence stars . Their resu l ts are 
presented in Tabled for S we give the mean abundance f rom lMorel et al. I (|2006l ) 



while the Ar abundance is taken from Lanz et al. ( 20081 ). When accounting for 



diffusion in the Sun, there is good agreement between the solar and hot star 
abundances for O, Ne and Mg in particular. While the situation has gr eatly 
improved compared with the solar abundances of lGrevesse &: Sauva 3 (Il998l ) and 
earlier compilations, there is still a tendency for the solar C, N and Fe abun- 
dances to be somewhat higher than in the B stars, contrary to expectations. It 
is unclear whether the solution can be found in the solar or B star analyses or, if 
a real difference indeed exists, perhaps due to infall of low-metallicity gas to the 
solar neighborhood. 

The abundances of nearby H 11 regions like the Orion nebula also provide con- 
straints on the solar c hemical compos i tion. In T able [5] we provide the values from 
Esteban et aD tooj ). lEsteban et all (j2005l ) and lGarcia-Roias fc EstebanI (|2007l ) 



interpolated to a Galacto-centric radius of 8 kpc. The observed scatter between 
individual Hii regions is now quite small and the inferred present-day Galactic 
abundance gradient well determined. It is worth remembering that there are 
multiple potential pitfalls in the analyses of H 11 regions, including corrections for 
unseen ionization stages and the effects of temperature fluctuations in the neb- 
ula. Furthermore, abundances derived from recombination (permitted) lines are 
often signiflcantly larger than the corresponding values from collisionally excited 
(forbidden) lines. Finally, dust condensation needs to be taken into account for 
some elements. For C and O this somew hat uncertain dust correction amounts to 
~ 0.1 dex in Orion ( Esteban et al.l 2004 ). Volatile elements like Ne and S are not 
expected to form dust but refractory elements like Mg and Fe are predominantly 
in the solid phase, so it is therefore not possible to deduce a meaningful total 
elemental abundance for them. Overall, the agreement between the Sun and H II 
regions is very satisfactory when accounting for diffusion and a minor chemical 
enrichment over the past 4.5 Gyr. 

Because of the generally lower temperatures encountered and the consequen- 
tially larger degree of dust condensation, the inte rstellar mediu m is a less reliable 
chemical composition gauge than Hii regions (Jenkins 20041 ) . Indeed, rather 
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than setting constraints on the solar abundances, the Sun is of ten used to in- 
fer th e amount in the sohd phase in the interstellar mediui n (e.g. ISofia &: Mevei 



2001). The local abundances of the volatile elements C ( Sofia et al. 2004). N 



( Jensen. Rachford &: SnoJ2007l ) and O (^ Cartledge et al ] l200d ] iJensen. Rachford &: Snow! 



2OO5I ) in the interstellar medium support our recommended solar values in gen- 
eral. With the old e r, hig h (log en = 8.83 — 8.93) solar ab undances preferred by 



Grevesse Sz Sauval ( 19981 ) and I Anders Grevesse ( 1989 ) it b ecomes very chal 



lenging to hide sufficient amounts of O in dust (Jensen 2007I ). The gas phase 
abundances of the refractory elements is only a small fraction of the total abun- 
dance in both the cold and warm interstellar medium, preventing them from 
setting any stringent constraints. 

Also listed in Table [5] are predict i ons fr om models of Galactic chemical evolution 
( Chiappini. Romano &: Matteucci 20031 ). Indeed, with these values it becomes 
difficult to reconcile the solar abundances with the results from OB stars and 
Hii regions for some elements like Fe. It is important to bear in mind though 
that the models have varying degrees of reliability for different elements. The 
evolutionary predictions for O, Ne and Mg should be the most accurate, as they 
are produced in hydrostatic burning and ejected by SNell; in fact for these ele- 
ments there is good agreement between the Sun and its local neighborhood. All 
other elements are either produced in explosive burning and therefore dependent 
on the adopted mass-cut in the supernova explosion, or have a non-negligible 
contribution from SNela or AGB stars, making theoretical estimates rather un- 
certain. One may conclude that the local Galactic volume has had a smaller 
star formation rate than assumed in these chemical evolution models. This in- 
ference is born out empirically by the o bserved age-metallicity relation, which 
iHolmberg. Nordstrom &: Andersen (I2OO8I ) estimated to be only 0.018 dex/Gyr for 
[Fe/H] in the solar neighborhood, i.e. only half of the predictions given in Table 
[5l An alternative explanation would be recent Galactic infall of more pristine 
gas. 



4.3 Helioseismology 

The situation is decidedly less encouraging when attempting to assess the so- 
lar abundances by means of helioseismology. Since the various solar oscillation 
p-modes penetrate to different depths, it is possible to thereby map the varia- 
tion of the sound speed in the solar interior, which can be com pared with the 
predicted values from sola r structure mode ls (see B asu fc Antial 2008, for a de- 
tailed review). With the Grevesse &: Sauvall (jl998l ) abundances, the predicted 
sound speed as a function of depth from standard solar models shows a very 
good agreeme nt with the helioseismic valu e s. W hen instead employing the re- 
vised values qf lAsplund. Grevesse &: Sauya l(l200,^). this is no longer the case (e.g . 



Bahcall et all l2005l. iBasu Antial l2004l. 1200 8'. 'Pelahave Pinsonneaultl [20oi 



Guzik. Watson Coxll2o'o5r iTurck-Chieze et al.ii2004i ). In this context it is per- 
haps fair to say that 'better is worse'. The sound speed discrepancy most clearly 
shows up immediately b elow the bottom of the c onvection zone, as illustrated 
in Fig. [HI Also with the Grevesse &: Sauval ( 19981 ) abundances there was a dis- 
agreement in this region but it is greatly aggravated by the further lowering of 
the solar O and Ne abundances and to a lesser extent by the reassessment of C, 
N and Fe. With the solar chemical composition recommended here, the situation 
is alleviated somewhat due to the slightly higher abundances of these elements 



28 



Asplund et al. 



and consequently larger opacity in the radiative interior but the deviation is still 
highly significant. 

It is not only the sound speed variation that is at variance when using the new 
solar abundances. The inferred depth of the convection zone is now too shallow: 
- Rrcz ^ 0.725 R(T) in stead of the helioseismic measurement 0.7133 ± 0.0005 R0 
(|Basu fc Antiall2004l l. Furthermore, the resulting He abundance when calibrat- 
ing solar interior models to achieve the correct solar luminosity and temperature 
at the solar age is similarly inconsistent with the value inferred from helioseis- 
mology ( see Sect. \3M : Ys « 0.238 compared with the observed 0.2485 ± 0.0034 
(jBasu Sz Antia 20041 ) . where Is is the present-day He mass fraction at the sur- 
face. The problem also persis ts in the solar core both as implied by low -degree 
p-modes ( Chaplin et al. 20071 ) and by gravity modes ( Garcia et al. 20071 ). Thus 
the problem exists both in the convective envelope, the radiative interior and in 
the core. 



Foll owin g the initial revisions of the solar C an d O abundances bv lAllende Prieto. Lambert 
(|200lh and kllende Prieto. Lambert fc Asplund (2002), many possible solutions 
have been put forward to explain the problem albeit to date with little success. 
The most obv i ous e xplanation would be that the opacities are underestimated. 
Bahcall et al. I (|2005l l estimated that a 10 — 20% increase of the OPAL opacities 
over a relatively substantial fraction of the solar i nterior (R = 0.4 — 0.7Rfri or 
equiva lently T = 2 — 5 • 10^ K) would be required; IChristensen-Dalsgaard et al. 



(I2OO9I) found even slightly more stringent r equir ements with the lAsplund. Grevesse &: Sauval 
(2005) solar abundances. Serenelli et al.l ( 20091 ) have investigated the differences 
with the solar chemical composition presented here and conclude that < 12% 
extra opacity is required to restore the previous agreement in sound speed. Lit- 
tle suggests, however, that current atomic physics calculations for s tellar interiors 



are rn issing such a substantial fraction of the opacity. For example, iBadnell et al 



showed that opacities from the Opacity Project are at most 3% larger 
than the corresponding OPAL values in the relevant temperature range. While 
unlikely to cure the problem by itself, it is perhaps wise to bear in mind that a 
similar situation for pulsating stars in fact existed in the 1980s, which found an 
unexpected resolution with the more comprehensive calculations by the OPAL 
and OP collaborations. 

There are other ways to compensate for the lower opacities resulting from the 
lower C, N, O and Ne abundances. One is that the radiative interior contains 
more metals tha n expected frora the p hotospheric abundances and standard mod- 
els for diffusion. lAsplund et al.l (j2004l ) first suggested that the diffusion efficiency 
may be seriously underestimated but the required elemental settlement is about 
twice the predicted value. Furthermore, while this may remedy the sound speed 
discrepancy, pro blems still remain in terms of the He abundance and depth of the 
conve ction zone ( Guzik. Watson &: Cox 2005, Montalban et al.l 2004 . Yang &: Bi 
2OO7I ). A hypothetical late accretion of met al-depleted gas during the forma- 



tion of the Sun runs into t he same trouble (jCastro. Vauclair &: RichardI 120071 : 
Guzik. Watson Coxl 20051 ) . Alternatively, the reduced O abundance could con- 



ceivabl y be compensated by a corresponding but larger i ncrea se in the Ne abun- 
dance (I Antia fc Basull2005l ^. iBahcall. Basu k Sereneilil (|2005l ^ estimated that a 
factor of three (0.5 dex) higher Ne content compared with the recommended val- 
ues bv lAsplund. Grevesse Sz Sauvall (j2005l ). together with minor adjustments of 
C, N, O and Ar within t heir associated unc e rtaint ies would restore the good 
agreement found with the Grevesse &: Sauval ( 19981 ) solar composition. There 
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was initial excitement when Drake &: Testa ( 20051 ) argued for exactly such an 



increase in Ne based on inferences from coronal spectroscopy of other stars but 
there is in fact little support for this idea, as explained in Sect. I3.9i The here 
recommended abundance of Ne as well as of C, N, O and Fe are ~ 0.05 dex 
higher than in Asplund. Grevesse &: Sauval ( 20051 ). This shift only partly allays 



the solar modelling discrepancy, as seen in Fig. [51 

Given that the problem exceeds the likely uncertainties in existing input data 
(e.g. opacities, diffusion, abundances), the only remaining option on the solar 
interior side seems to be to invoke some physical process not yet accounted for 
in standard solar models. It is perhaps suggest ive that the largest deviat i on oc - 
curs immediately below the convection zone. lArnett. Meakin &: Youn3 (|2005l ) 



have proposed that internal gravity waves generated by the convective motion 
at the base of the convection zone propagate inwards and deposit their energy 
in the radiative zone. The ensuing structure changes are similar to an opacity 
enhancement, thus going in the right direction. I n addition, the g r avity waves 
induce additional mixing, which would also help. Young &: Arnett ( 20051 ) have 



investigated these processes using multi-dimensional hydro dynamical simulations 
of stellar convection zones and found qualitative improvements in many stellar 
regimes. Due to the computationally very challenging nature of the hydrodynam- 
ics, the modelling is unfortunately not yet able to quantitatively predict whether 
the effects are sufficient to explain the helioseismic signature. 

Thus, most suggestions to resolve the solar modelling discrepancy have already 
been ruled out or are considered improbable. The few possibilities still stand- 
ing do so largely because no detailed modelling has as yet been forthcoming to 
properly put the proposals to test. It is possible though that a combination of 
several factors could be the explanation, but such fine-tuning would seem rather 
contrived. Should no likely solution be found, it would suggest that the photo- 
spheric abundances presented herein are wrong. Given the discussion in Sect. [2] 
and Sect. [Sj it is not easy to understand where such a serious error may lie but 
the possibility should obviously not be discounted. 

4.4 Solar neutrinos 

Neutrinos provide another means by which to infer the solar chemical compo- 
sition, though current experiments are not yet sufficiently sensitive to return 
decisive results in this respect. The metallicity adopted in a solar model mod- 
ifies the conditions in the core, and thus the neutrino fluxes produced by the 
pp-chain and CNO-cycle. Whilst the neutrino fluxes cannot be directly asso- 
ciated with abundances of individual element s, the relevant neutrino c hannels 
depend most sensitively on C, N, O, Si and Fe (|Bahcall fc Serenellil[2005l ). Using 



neutrino data from the S uperK amiokande I-I-II, SNO and Borexino experiments, 
Pena-Garav &: Serenelli (|2008l l have derived the neutrino oscillations and fluxes 



expected from solar models c onstr ucted with the Grevesse &: Sauval ( 19981 ) and 



Asplund. Grevesse &: Sauva 3 (|2005l ) compositions. Intriguingly, the existing data 



imply ^Be and ®B neutrino fluxes that fall between the predictions for the two 
models; we suspect that the solar neutrinos would favour the new abundances 
presented herein. A more direct test of the heavily debated solar C, N and 
O abundances would come from measuring the ^^N and ^^O /3-decay neutrinos 
that follow proton capture on ^^C and ^'^N, respectively. Meaningful constraints 
should come from Borexino within the next couple of years, and the proposed 
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SNO+ experiment should weigh in with a definite answer sometime after 2011 
(jHaxton fc Serenellil[2008l l. 



5 CONCLUDING REMARKS 
5.1 Summary Points 

We have critically examined all ingredients for the determination of the solar pho- 
tospheric chemical composition, from the atmospheric and line formation mod- 
elling, required atomic data and selection of lines to the final estimate of the 
elemental abundances and associated errors (Sect. [3]). A notable feature of this 
analysis is that for essentially every element, the effects of 3D hydrodynamical 
model atmospheres have now been accounted for. The 3D solar model employed 
here outperforms any ID model, theoretical or semi-empirical, in a raft of obser- 
vational tests. The values we recommend for the most abundant metals C, N, O, 
Ne and Fe are significantl y smaller (by ^ 0.2 dex) th an those advocated in the 
widely-used compilation of Anders &: Grevesse ( 19891 ). Most of these differences 



can be attributed to the improvements over the past two decades in modelling 
the solar atmosphere, allowing for departures from LTE in the line formation, 
properly treating blends and employing improved atomic and molecular data 
(Sect. [2]). Amongst other things, these new, lower abundances are supported by 
the excellent agreement between different abundance indicators for C, N and O: 
low-excitation forbidden atomic lines, high-excitation permitted atomic lines and 
various molecular features. 

A comparison of the photospheric abundances and the most primitive class of 
meteorites, the CI chondrites, reveals a very good agreement for nearly all non- 
volatile elements (Sect. 14. ip . Since the volatile elements, including H, have been 
depleted in the chondrites, the meteoritic abundance scale has been placed on the 
same absolute level as the photospheric abundances by enforcing equivalence of 
the Si abundances. The remaining few elements (Rb, Rh, Pd, Ag, Hf, Pb) that 
suggest a significant difference between meteorites and the Sun can probably 
be attributed to poor solar determinations, likely due to incorrect ^/-values, 
unidentified blends or unaccounted non-LTE effects. 

In order to make a meaningful comparison between the Sun and the abundances 
measured in the solar neighborhood, one must first correct for the effects of 
diffusion in the Sun and chemical enrichment in the Galaxy over the past 4.56 Gyr. 
All elements are predicted to have settled to the solar interior by essentially the 
same amount, meaning that the bulk, or equivalently the proto-solar, composition 
is higher by 0.04 dex than listed in Table [T] for all elements heavier than He 
(Sect. 13. lip . Our recommended solar abundances agree very well in general with 
those measured in nearby OB stars, Hii regions like Orion and the local ISM 
(Sect. 14. 2[) . We note that care should be exercised when comparing our solar 
results with corresponding abundances for late-type stars derived with classical 
ID model atmospheres; ideally a differential study of the stars and the Sun should 
be performed and using the solar abundances presented in Tabled] to place them 
on an absolute scale. 

In spite of the many apparent successes of the solar chemical composition we 
recommend here, the outlook is not entirely rosy. With the downward revision 



of the overall metal content presented here compared to the iGrevesse Sz Sauval 



( 19981 ) abundances, the previous excellent agreement between predicted sound 
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speeds from standard solar models and those inferred from helioseismology is 
ruined (Sect. 14. Sp . To date there has been no fully convincing solution put 
forward. The dis cordance has been alleviated some what relative to the recom- 
mended values in Asplund. Grevesse Sauval ( 20051 ) but it nevertheless remains 



a significant discrepancy in urgent need of resolution. 
5.2 Future Issues 

While major progress has been made in recent years in the inferred solar chem- 
ical composition, there is no question that a great deal of work remains before 
we can confidently state that each solar elemental abundance is known to bet- 
ter than 0.01 dex (2%). For many, if not most elements, the achieved accuracy 
today is more like 0.05 dex or indeed significantly worse in many cases. A con- 
certed effort on multiple fronts is required to substantially decrease the remaining 
uncertainties. 

We do not consider it likely that subsequent generations of 3D hydro dynam- 
ical model atmospheres will dramatically (i.e. at the > 0.1 dex level) amend 
the inferred solar photospheric chemical composition, given how successful the 
most recent models are in reproducing key observational diagnostics. It is fair 
to say that the mean temperature structure of the quiet solar granulation, the 
velocity distribution and the atmospheric inhomogeneities are now rather well 
known, at least in the formation regions of the weaker lines largely employed by 
the abundance studies described herein. Higher atmospheric layers, to which the 
strongest atomic and molecular lines are more sensitive, are naturally less con- 
strained observationally; on the modelling side this requires improved radiative 
transfer treatment and allowance of departures from LTE for key elements. For 
completeness, one should include the effects of magnetic fields in the 3D solar 
atmospheric models, bu t we do not expect that these wil l noticeably influence 
the derived abundances ( Nordlund. Stein &: Asplundl 20091 ). 



The greatest uncertainties afflicting solar abundance work today instead stem 
from poor atomic and molecular data, and incomplete modelling of departures 
from LTE. Few spectral lines have their transition probabilities determined to bet- 
ter than a few %. As an example here it suffices to mention the case of Mg with 
its relatively simple atomic structure: the majority of the Mgi lines employed 
here and elsewhere for solar abundance purposes still have quoted uncertainties 
of 25 — 50% according to the NIST database, which is clearly unacceptable. For- 
tunately there has been some laudable pro gress on this score in recent years 



prog 

both in terms of experime i ital efforts fe.g. iBlackwell- Whitehead Sz Bergemann 
2007 . Johanss on et al. '2003','Lawler et al."2009'), and theoretical calculations (e.g. 



Badnell et al.i, 2005, . .Kurucz . 1992.. Plez, 1998) as discussed in Sect. 12.21 A major 
problem, however, is the dearth of atomic physicists devoted to satisfying the 
astronomical need for various types of input data; the astronomical community 
ought to better appreciate and encourage the relatively few individuals carrying 
out this important work, given its tremendous benefit for solar and stellar as- 
trophysics especially. On the wish-list are improved 5/- values for Mgl, Mgll, 
Sin, Call, the once-ionized species of the Fe-peak elements, Cui and Srii. Also, 
atomic lines in the IR have great potential given the accurate determination of 
the continuum and few blends but more often than not accurate transition prob- 
abilities are lacking. 

Equally important is the consideration of non-LTE effects on line formation. 
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As discussed by Asplund ( 20051 ). even though departures from LTE are expected 
to be relatively serene in the Sun, one should always expect non-LTE abundance 
corrections at around the 0.05 dex level; indeed in many cases the effects will be 
significantly more severe. Still today, most elements have not been exposed to a 
detailed non-LTE study even in the Sun. Even rarer are solar abundance works 
going beyond the ID framework and employing 3D no n-LTE calculations (e.g. 
Asplund. Carlsson fc Botnenll2003l . I Asplund et"alll2004l ^. In part, the scarcity of 



non-LTE studies is related to a lack of adequate atomic data. Compared with the 
LTE case, vastly more input data are required in non-LTE calculations than are 
presently available for most elements. Significant headway will only become pos- 
sible with the advent of detailed photo-ionization and collisional cross-sections; 
for the latter, estimates of the impact of both electron and hydrogen collisions, as 
well as other processes like charge transfer, are urgently needed. There is also an 
issue of manpower; careful non-LTE studies are very time-consuming and there 
are not enough specialists working in the field to satisfy the demand. Finally, we 
note that non-LTE line formation should be explored in 3D model atmospheres, 
since the the presence of atmospheric inhomogeneities is generally expected to 
amplify any existing non-LTE effects seen in ID. 

A section on future issues in the solar chemical composition would not be com- 
plete without stressing the importance of finally resolving the solar modelling 
problem uncovered by helioseismology (Sect. 14. 3p . Many hypotheses have been 
put forward to explain the discrepancy, but no clear solution has been forthcom- 
ing. Until an explanation has been found, a mutual cloud of doubt will always 
linger over the solar chemical composition recommended here and current models 
of the solar interior. Whether the problem will be traced in the end to shortcom- 
ings in the atmospheric or interior modelling remains to be seen, but whatever 
the resolution it will place solar and stellar astrophysics - and by extension most 
of astronomy - on a much firmer footing. 
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Table 1: Element abundances in the present-day solar photosphere. 
Also given are the corresp onding values for CI carbonaceous chondrites 



( Lodders. Palme &: Gaill 120 03). Indirect photospheric estimates have been used 



for the noble 


gases (Sect. 13.91 














Elem. 


Photosphere 


Meteorites 




Elem. 


Photosphere 


Meteorites 
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12.00 


8.22 ±0.04 


44 


Ru 


1.75 ±0.08 


1.76 ±0.03 
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He 


[10.93 ±0.01] 


1.29 


45 


Rh 


0.91 ±0.10 


1.06 ± 0.04 
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Li 


1.05 ±0.10 


3.26 ±0.05 


46 


Pd 


1.57 ±0.10 


1.65 ±0.02 
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Be 


1.38 ±0.09 


1.30 ±0.03 


47 


Ag 


0.94 ±0.10 


1.20 ±0.02 
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2.70 ±0.20 


2.79 ±0.04 


48 


Cd 




1.71 ±0.03 
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8.43 ±0.05 


7.39 ±0.04 


49 


In 


0.80 ±0.20 


0.76 ± 0.03 
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7.83 ±0.05 


6.26 ±0.06 


50 


Sn 


2.04 ±0.10 


2.07 ±0.06 
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Te 
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Ne 
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11 


Na 


6.24 ± 0.04 


6.27 ±0.02 


54 


Xe 


[2.24 ± 0.06] 


-1.95 


12 


Mg 


7.60 ± 0.04 


7.53 ±0.01 


55 


Cs 




1.08 ±0.02 


13 


Al 


6.45 ± 0.03 


6.43 ±0.01 


56 


Ba 


2.18 ±0.09 


2.18 ±0.03 


14 


Si 


7.51 ±0.03 


7.51 ±0.01 


57 


La 


1.10 ±0.04 


1.17 ±0.02 


15 


P 


5.41 ± 0.03 


5.43 ±0.04 


58 


Ce 


1.58 ±0.04 


1.58 ±0.02 


16 


S 


7.12 ± 0.03 


7.15 ±0.02 


59 


Pr 


0.72 ±0.04 


0.76 ± 0.03 


17 


CI 


5.50 ±0.30 


5.23 ±0.06 


60 


Nd 


1.42 ±0.04 


1.45 ±0.02 


18 


Ar 


[6.40 ±0.13] 


-0.50 


62 


Sm 


0.96 ±0.04 


0.94 ±0.02 


19 


K 


5.03 ±0.09 


5.08 ±0.02 


63 


Eu 


0.52 ±0.04 


0.51 ±0.02 


20 


Ca 


6.34 ± 0.04 


6.29 ±0.02 


64 


Gd 


1.07 ±0.04 


1.05 ±0.02 


21 


Sc 


3.15 ±0.04 


3.05 ±0.02 


65 


Tb 


0.30 ±0.10 


0.32 ±0.03 


22 


Ti 


4.95 ± 0.05 


4.91 ±0.03 


66 


Dy 


1.10 ±0.04 


1.13 ±0.02 


23 


V 


3.93 ± 0.08 


3.96 ±0.02 


67 


Ho 


0.48 ±0.11 


0.47 ±0.03 


24 


Cr 


5.64 ± 0.04 


5.64 ±0.01 


68 


Er 


0.92 ±0.05 


0.92 ±0.02 


25 


Mn 


5.43 ± 0.05 


5.48 ±0.01 


69 


Tm 


0.10 ±0.04 


0.12 ±0.03 


26 


Fe 


7.50 ± 0.04 


7.45 ± 0.01 


70 


Yb 


0.84 ±0.11 


0.92 ± 0.02 


27 


Co 


4.99 ± 0.07 


4.87 ±0.01 


71 


Lu 


0.10 ±0.09 


0.09 ± 0.02 


28 


Ni 


6.22 ±0.04 


6.20 ±0.01 


72 


Hf 


0.85 ±0.04 


0.71 ± 0.02 


29 


Cu 


4.19 ±0.04 


4.25 ± 0.04 


73 


Ta 




-0.12 ±0.04 


30 


Zn 


4.56 ±0.05 


4.63 ± 0.04 


74 


W 


0.85 ±0.12 


0.65 ± 0.04 


31 


Ga 


3.04 ± 0.09 


3.08 ±0.02 


75 


Re 




0.26 ± 0.04 


32 


Ge 


3.65 ±0.10 


3.58 ±0.04 


76 


Os 


1.40 ±0.08 


1.35 ±0.03 


33 


As 




2.30 ± 0.04 


77 


Ir 


1.38 ±0.07 


1.32 ±0.02 


34 


Se 




3.34 ± 0.03 


78 


Pt 




1.62 ±0.03 


35 


Br 




2.54 ±0.06 


79 


Au 


0.92 ±0.10 


0.80 ± 0.04 


36 


Kr 


[3.25 ± 0.06] 


-2.27 


80 


Hg 




1.17 ±0.08 


37 


Rb 


2.52 ±0.10 


2.36 ±0.03 


81 


Tl 


0.90 ±0.20 


0.77 ±0.03 


38 


Sr 


2.87 ±0.07 


2.88 ±0.03 


82 


Pb 


1.75 ±0.10 


2.04 ±0.03 


39 


Y 


2.21 ±0.05 


2.17 ±0.04 


83 


Bi 




0.65 ± 0.04 


40 


Zr 


2.58 ±0.04 


2.53 ±0.04 


90 


Th 


0.02 ±0.10 


0.06 ± 0.03 


41 


Nb 


1.46 ± 0.04 


1.41 ±0.04 


92 


U 




-0.54 ± 0.03 


42 


Mo 


1.88 ±0.08 


1.94 ±0.04 
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Table 2: The derived photospheric solar C, N and O abundances as indicated 
by various atomic and molecular species a nd a 3D hydr odynamical model of the 
solar atmosphere ( Asplund et alJ 2009al 3. Sauval et al. 2009. ) . The correspond- 
ing results obtained with the temporally and spatially (over equal optical depth 
surfaces) averaged 3D m odel (denoted < 3D>), the ID theor etical MARCS model 
Gustafsson et al. 20081 ) and the ID semi-empirical model of Holweger fc Miiller 



I974I ) are also given for comparison purposes. The errors given here are only the 



standard deviation (cr) for each type of transition. For the atomic transitions, 
departures from LTE have been considered. 



lines 






logec,N,o 








3D 


< 3D> 


HM 


MARCS 


[Ci] 




8.41 


8.40 


8.41 


8.38 


Ci 




8.42 ±0.05 


8.47 ± 0.04 


8.45 ±0.04 


8.39 ± 0.04 


CH Av = 


1 


8.44 ± 0.04 


8.44 ± 0.04 


8.53 ±0.04 


8.44 ± 0.04 


CH A-X 




8.43 ± 0.03 


8.42 ± 0.03 


8.51 ±0.03 


8.40 ± 0.03 


C2 Swan 




8.46 ± 0.03 


8.46 ± 0.03 


8.51 ±0.03 


8.46 ± 0.03 


Ni 




7.78 ± 0.04 


7.89 ± 0.04 


7.88 ±0.04 


7.78 ± 0.04 


NH Av = 





7.83 ±0.03 


7.94 ± 0.02 


8.02 ±0.02 


7.97 ± 0.02 


NH Av = 


1 


7.88 ± 0.03 


7.91 ±0.03 


8.01 ±0.03 


7.91 ±0.03 


[Oi] 




8.70 ± 0.05 


8.70 ± 0.05 


8.73 ±0.05 


8.69 ± 0.05 


Oi 




8.69 ± 0.05 


8.73 ± 0.05 


8.69 ±0.05 


8.62 ± 0.05 


OH Av = 





8.69 ±0.03 


8.75 ±0.03 


8.83 ±0.03 


8.78 ±0.03 


OH Av = 


1 


8.69 ±0.03 


8.74 ± 0.03 


8.86 ±0.03 


8.75 ± 0.03 
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Table 3: Representative isotopic abu ndance fractions in the s olar system. Most 
of the isotopic values are taken from Rosman &: Tavlor ( 19981 ) with updates for 
some elements, as discussed in Sect. 13.101 
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95 
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Ne 20 
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47 
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74 


36.28 




96 


16.672 
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21 


0.2228 




48 


73.72 


76 


7.61 




97 


9.599 




12032.58 




22 


6.8341 




49 


5.41 








98 


24.391 




1224.63 










50 


5.18 


As 75 


100.0 




1009.824 




1245.79 


Na 23 


100.0 




























V 


50 


0.250 


Se 74 


0.89 


Ru 96 


5.54 


Sb 


12157.21 


Mg 


24 


78.99 




51 


99.750 


76 


9.37 




98 


1.87 




12342.79 




25 


10.00 








77 


7.63 




99 


12.76 








26 


11.01 


Cr 


50 


4.345 


78 


23.77 




10012.60 


Te 


1200.09 










52 


83.789 


80 


49.61 




10117.06 




1222.55 


Al 


27 


100.0 




53 


9.501 


82 


8.73 




10231.55 




1230.89 










54 


2.365 








10418.62 




1244.74 


Si 


28 


92.2297 








Br 79 


50.69 










125 7.07 




29 


4.6832 


Mn55 


100.0 


81 


49.31 


Rh 103100.0 




12618.84 




30 


3.0872 




















12831.74 








Fe 


54 


5.845 


Kr 78 


0.362 


Pd 1021.02 




13034.08 


P 


31 


100.0 




56 


91.754 


80 


2.326 




10411.14 
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14913.82 


16722.869 




18628.43 




2046.87 


Cs 


133100.0 


1507.38 


16826.978 














15226.75 


17014.910 


Re 18535.662 


Tl 


20329.524 


Ba 1300.106 


15422.75 






18764.338 




20570.476 




1320.101 




Tml69 100.0 












1342.417 


Eu 15147.81 




Os 1840.020 


Pb 2041.997 




1356.592 


15352.19 


Yb 1680.12 




1861.598 




20618.582 




1367.854 




1702.98 




1871.271 




20720.563 




13711.232 


Gd 1520.20 


17114.09 




18813.337 




20858.858 




13871.698 


1542.18 


17221.69 




18916.261 










15514.80 


17316.10 




19026.444 


Bi 


209100.0 


La 


1380.091 


15620.47 


17432.03 




19241.070 








13999.909 


15715.65 


17613.00 






Th 232 100.0 






15824.84 




Ir 


19137.3 






Ce 


1360.185 


16021.86 


Lu 17597.1795 




19362.7 


U 


2340.002 




1380.251 




1762.8205 








23524.286 




14088.450 


Tb 159100.0 




Pt 


1900.014 




23875.712 




14211.114 




Hf 1740.162 




1920.782 










Dy 1560.056 


1765.206 




19432.967 






Pr 


141 100.0 


1580.095 


17718.606 




19533.832 
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Table 4: The mass fractions of hydrogen (X), hehum (Y) and metals (Z) for a 
number of widely-used compilations of the solar chemical composition. 



Source 




X 




Y 




z 


z/x 


Present-day photosphere: 


















Anders k Grevesse fl989)'' 





7314 





2485 





0201 





0274 


Grevesse & Noels ri993)^ 





7336 





2485 





0179 





0244 


Grevesse k Sauval ^1998) 





7345 





2485 





0169 





0231 


Lodders f 20031 





7491 





2377 





0133 





0177 


AsDlund. Grevesse k Sauval (2005) 





7392 





2485 





0122 





0165 


Lodders. Palme k Gail (2009) 





7390 





2469 





0141 





0191 


Present work 





7381 





2485 





0134 





0181 



Proto-solar: 



Anders k Grevesse fl989) 





7096 





2691 





0213 





0301 


Grevesse k Noels fl993) 





7112 





2697 





0190 





0268 


Grevesse k Sauval ^1998) 





7120 





2701 





0180 





0253 


Lodders (2003) 





7111 





2741 





0149 





0210 


AsDlund, Grevesse k Sauval (2005) 





7166 





2704 





0130 





0181 


Lodders. Palme k Gail f2009) 





7112 





2735 





0153 





0215 


Present work 





7154 





2703 





0142 





0199 



^ The He abundances given in Anders k Grevesse (|l989l ) and lGrevesse k Noelsl 
( I993I ) have here been replaced with the current best estimate from helioseismol- 
ogy (Sect. [3:911 . 



Table 5: Comparison of t h e pro to-solar abundances from the present 
work and Grevesse k Sauval ( 19981 ) with those in nearby B stars and 
Hii regions. The solar values given h ere include the effects of diffusion 
( Turcotte k Wimmer-Schweingrubeill2002l ) as discussed in Sect. 13.111 The Hii 
numbers include the estimated elemental fractions tied up in dust; the dust cor- 
rections for Mg, Si and Fe are very large and thus too uncertain to provide 
meaningful values here. Also given in the last column is the predicted Galactic 
chemical enrichment (GCE) over the past 4.56 Gyr. 

Elem. Sun^ Sun^ B stars'^ H ii^^ GCE^ 



He 


10.98 ±0.01 


10.98 ±0.01 


10.98 ± 0.02 


10.96 ±0.01 





01 


C 


8.56 ± 0.06 


8.47 ±0.05 


8.35 ±0.03 


8.66 ± 0.06 





06 


N 


7.96 ± 0.06 


7.87 ±0.05 


7.76 ± 0.05 


7.85 ± 0.06 





08 





8.87 ±0.06 


8.73 ±0.05 


8.76 ± 0.03 


8.80 ± 0.04 





04 


Ne 


8.12 ±0.06 


7.97 ±0.10 


8.08 ± 0.03 


8.00 ± 0.08 





04 


Mg 


7.62 ± 0.05 


7.64 ±0.04 


7.56 ± 0.05 







04 


Si 


7.59 ± 0.05 


7.55 ±0.04 


7.50 ± 0.02 







08 


S 


7.37 ±0.11 


7.16 ±0.03 


7.21 ±0.13 


7.30 ± 0.04 





09 


Ar 


6.44 ± 0.06 


6.44 ±0.13 


6.66 ± 0.06 


6.62 ± 0.06 






Fe 


7.55 ± 0.05 


7.54 ±0.04 


7.44 ± 0.04 







14 



Grevesse k Sauval ( 19981) ^ Present work ^ Przvbilla. Nieva k Butlei] 


(2008). 


Morel etal. (2006). Lanz et al. (2008) ^ Esteban et al. 


(2005. I2OO4I). 


Garcia-Rojas k Esteban ( 


2OO7I) ^ IChiappini, Romano k Matteucci 


(20031). 
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Fi gure 1: The mean temp erature structure of the 3D hydrodynamical model 
of IXrampedach et al. ( 20091 ) is shown as a function of optical depth at 500 nm 
(blue solid line). The blue dashed lines correspond to the spatial and tem- 
poral rms variations of the 3D model, while th e red and green curves denote 
the ID semi-empirical iHolweger fc Miiller ( 19741 ) and the ID theoretical marcs 
( Gustafsson et al. 20081 ) model atmospheres, respectively. 
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Figure 2: Comparison of the predicted continuum center-to- limb variation as a 
function of wavelengt h for different solar model atmospheres against observations 
(iNeckel Labslll994l ). The results for five different viewing angles fi are shown 
from near disk center {fj, = 0.9) to close to the limb (/x = 0.1). All intensities are 
normal ized to the correspondin g disk-center intensitie s. The 3D hvdrodynamica l 
model ( Trampedach et 311120091 ) outperforms even the lHolweger k Miilleil (Il974l) 



semi- empirical model, which was designed to satisfy this diagnostic (iPereira et al 
20091 ). As for all ID theoretical model atmospheres, the marcs model has a too 



steep temperature gradient, which manifests itself in a poor agreement with the 
center-to-limb variation. 
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Figure 3: The predicted spectral hne profile of a typical Fei line from the 3D 
hydrodynamical solar model (red solid line) compared with the observations (blue 
rhombs). The agreement is clearly very satisfactory, which is the result of the 
Doppler shifts arising from the self-consistently computed convective motions that 
broaden, shift and skew the theoretical profile. For c omparison purposes also the 
predicted profile from a ID model atmosphere (here Holweger &: Miiller 19741 ) is 
shown; the ID profile has been computed with a microturbulence of 1 kms~^ and 
a tuned macroturbulence to obtain the right overall linewidth. Note that even 
with these two free parameters the ID profile can neither predict the shift nor 
the asymmetry of the line. 
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20 40 60 80 

Atomic number 

Figure 4: The present-day solar photospheric elemental abundances as a function 
of atomic number. As throughout this article, the logarithmic abundance of 
hydrogen is defined to be log en = 12.0. 
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Figure 5: Comparison of the observed (open circles with error bars) and predicted 
center-to-hmb variation of the equivalent width of the Ol 777.19 nm triplet line 
as a function of viewing angle (/x = 1 corresponds to disk center and to the 
limb). The theoretical curves are the results of full 3D non-LTE line forma- 
tion calculations with different adopted efficiency for the poorly known inelas- 
tic H collisions; the O abundances have been normalized to fit the disk-center 
line strengths in each case. The curves ar e labelled with the scaling factor 
used to multiply the classical Drawin (jl968l ) formula. The solar observations 
thus seem to suggest that this recip e is a rather accurate representation for O I 
( Pereira. Asplund &: Kiselman 20091 ). 
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Figure 6: Upper panel: The derived Fei (blue circles) and Fell (red circles) 
abundances ( Scott et alJ 2009bl ) as a function of excitation potential of the lower 
level when using a 3D hydrodynamical solar model atmosphere. Lowe r panel- 
Same as above but with the ID semi-empirical iHolweger &: Mullerl \W74 ) model. 
Note the offset between Fe i and Fe 11 and the trend with excitation potential for 
the Fe l lines, neither of which are present when employing the 3D model. 
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Figure 7: Difference between the logarithmic abundances determined from the 
solar photosphere and the CI carbonaceous chondrites as a function of atomic 
number. With a few exceptions the agreement is excellent. Note that due to 
depletion in the Sun and meteorites, the data points for Li, C, N and the noble 
gases fall outside the range of the figure. 
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Figure 8: The difFeren ces between the helio seismic and predicted sound speeds 
as a function of depth ( Serenehi et al. 20091 ). The standard solar models shown 



here only differ in the assumed ch emical compositions: [ Grevesse & Sauva 
(black line, here denoted GS98), lAsplund. Grevesse fc Sauval (20q1) (red line, 
AGS05) and the present work (blue line, AGSS09). Each model has indepen- 
dently been calibrated to achieve the correct solar luminosity, temperature and 
age. The base of the convection zone is at = 0.71 R0, which is also where the 
discrepancy starts in earnest in all three cases. 



